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Abstract 



We study the fluxes of neutrinos from annihilations of dark matter parti- 
cles in the Sun and the Earth. We give the spectra of all neutrino flavors 
for the main known annihilation channels: z/z>, 66, rf, cc, light quarks, 
ZZ, W^W^ . We present the appropriate formalism for computing the 
combined effect of oscillations, absorptions, i/T--regeneration. Total rates 
are modified by an 0{Q.l 10) factor, comparable to astrophysical un- 
certainties, that instead negligibly affect the spectra. We then calculate 
different signal topologies in neutrino telescopes: through-going muons, 
contained muons, showers, and study their capabilities to discriminate 
a dark matter signal from backgrounds. We finally discuss how mea- 
suring the neutrino spectra can allow to reconstruct the fundamental 
properties of the dark matter: its mass and its annihilation branching 
ratios. 



1 Introduction 



The most appealing scenario to explain the observed Dark Matter (DM) abundance ^Ibm ~ 0.3 
consists in postulating that DM arises as the thermal relic of a new stable neutral particle with 
mass mDM- Assuming it has weak couplings ~ 1, the right ^dm is obtained for ttidm ~ 
(TMpi)^/^ ~ TeV, where T ~ 3 K is the present temperature of the universe, and Mpi ~ 10^^ GeV 
is the Planck mass [T] . One motivated DM candidate is the lightest neutralino in supersymmetric 
extensions of the Standard Model with conserved matter parity, that for independent reasons 
is expected to have a mass around the electroweak scale [2]. Many other DM candidates have 
been proposed: we will generically have in mind a DM particle heavier than few tens of GeV, 
keeping the concrete connection to the neutralino as a guideline. This scenario seems testable 
by DM search and by collider experiments: one would like to see a positive signal in both kind 
of experiments and to check if the same particle is responsible for both signals. As emphasized 
in [3] this is an important but difficult goal. 

A huge effort is currently put in experiments that hope to discover DM either directly (through 
the interaction of DM particles with the detector) or indirectly (through the detection of secondary 
products of DM annihilations). Among the indirect methods, a promising signal consists in 
neutrinos with energy Ei, < rriDM produced by annihilations of DM particles accumulated in the 
core of the Earth and of the Sun [H [5], detected by large neutrino detectors. We will refer 
to them as 'DMz/'. 1MB [6], Kamiokande [7], Baksan [8], Macro [9], Super-Kamiokande [10], 
AMANDA [11] and BAIKAL [12] already obtained constraints on DMu fluxes, while experiments 
that are under construction, like ANTARES [13] and ICECUBE [2], or that are planned, like 
NEMO [T3], NESTOR [H] and a Mton-scale water Cerenkov detector [T7], will offer improved 
sensitivity. 

We compute the spectra of neutrinos of all flavors generated by DM annihilations in the Earth 
and in the Sun. 

Today, before a discovery, this can be used to convert experimental data into more reliable 
constraints on model parameter space and helps in identifying more relevant features of the DMi^ 
signal searched for. For instance, we include in the analysis all main annihilation channels, we 
address the effect of neutrino oscillations and interactions with matter and we point out more 
experimental observables that those usually considered. 

After a discovery the situation will be analogous to the solar neutrino anomaly: a natural 
source of neutrinos carries information about fundamental parameters and we must find realistic 
observables that allow to extract it. As in that case, also in the DM case the total u rate is the 
crucial parameter for discovery but is plagued by a sizable (9(10) astrophysical uncertainty. How 
can we then reconstruct the properties of the DM? 

Astrophysical uncertainties negligibly affect the ratios between different neutrino flavors and 
the neutrino energy spectra (as well as the closely related angular distributions [IS])- They depend 
on the DM mass mDM and on the branching ratios of the channels into which DM particles may 
annihilate: uu, bb, tt, t^t~ , W'^W" , ZZ... In order to extract these fundamental parameters from 
future data one needs to precisely compute DMi/ spectra taking into account the astrophysical 
environment, where several processes are important. 

In section |2] we motivate our phenomenological procedure and compute the fluxes of electron. 
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Figure 1: The left plot illustrates how oscillations and CC absorption separately affect a flux of 
neutrinos produced in the center of the Sun. The right plot shows the oscillation probabilities from 
the center of the Earth. The continuous line applies to v for 613 = and to v for any allowed ^13, 
since matter effects suppress their mixing. The dotted line applies to v for 6'i3 = 0.1 rad. The 
average over the production point has been performed as appropriate for ttidm = 100 GeV. It is 
responsible for the damping effect visible at Ey < 10 GeV. 



muon and tau neutrinos at production point: the different density of the Earth and solar core 
affects energy loss of particles that decay producing neutrinos. 

In section [3] we compute how propagation from the center of the Earth and of the Sun affects 
the flavor and energy spectra. At production, the neutrino flavor ratios from the DM annihilations 
are simply given by: 

z/e : z/g : i/^ : z/^ : i/^ : z/^ = 1 : 1 : 1 : 1 : r : r . 

In the Earth, the main effect is due to oscillations with 'atmospheric' frequency: the neutrino 
oscillation length 

Aatm = inE,/\AmlJ ^ 10^ km(E,/100 GeV) 

is comparable to the Earth radius = 6371 km if E^, < 100 GeV. In the Sun, also the size of the 
production region of DMz/ is of the same order. Furthermore, in the Sun at E^ > 10 GeV neutrino 
interactions start to be significant and solar oscillations cease to be adiabatically MSW-enhanced, 
as illustrated in fig.[T^. Interactions manifest in several ways: absorption, re-injection of neutrinos 
of lower energy (as produced by NC scatterings and Ur CC scatterings), breaking of coherence 
among different flavors. These effects operate at the same time and with comparable importance: 
while previous works addressed the issues separately [ISl [20], the density-matrix formalism of 
section [3] allows to take into account their combined action. 
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In section |4] ([5]) we give the resulting energy spectra of DM neutrinos of all flavors from the 
Earth (Sun). We consider the standard through-going muon signal and point out that other 
classes of events can be studied in realistic detectors and have interesting features from the point 
of view of discriminating a DMu signal from the atmospheric background and of reconstructing 
DM properties. This latter point is discussed in section |6l 



2 Neutrino production 



A flux of neutrinos is produced inside the Earth or the Sun as a consequence of annihilation 
of dark matter particles which have been gravitationally captured inside these celestial bodies 
The differential neutrino flux is: 



dE^ And? ^"^^^'dE 



where / runs over the different final states of the DM annihilations with branching ratios BR/, 
d is the distance of the neutrino source from the detector (either the Sun-Earth distance tse or 
the Earth radius R^) and where the annihilation rate Fann depends on the rate Fcapt of captured 
particles by the well known relation: 

rann = ^tanh^ (Vta) (2) 

where to = 4.5 Gyr is the age of the Earth and of the Sun and ta denotes a time-scale for the 
competing processes of capture and annihilation, and it is proportional to the DM annihilation 
cross sections (for explicit formulae see [211 l23l El])- For the present discussion we just remind 
that the capture rate F^apt depends linearly on the DM/nucleus scattering cross section and on 
the local dark matter density Pdm: 

^ capt ''"scattering Pdm (3) 

Eq. ([2]) shows that the two competing processes of capture and annihilation may eventually reach 
an equilibrium situation when the time scale ta is much smaller than the age of the body. While 
this is usually the case for the Sun, it does not always occur for the Earth, since in this case 
the gravitational potential, which is responsible for the capture, is much smaller. Equilibrium is 
fulfilled only for large elastic scattering cross sections. 



2.1 Observables with and without astrophysical uncertainties 

From the previous equations we see that the neutrino signal shares both astrophysical and particle 
physics uncertainties. However, the shape of neutrino spectra are virtually free from the astro- 
physical ones, even in presence of oscillations - as we shall discuss below - and therefore they 
can be potentially used to study the fundamental DM parameters, like its mass and annihilation 
channels. This topic will be addressed in section El 

A quantity which suffers from sizable astrophysical uncertainties is the total DMz/ flux, mainly 
due to the poor knowledge of the local DM density pdm- The experimental indetermination on 
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this parameter is still large. Detailed analyses, performed assuming different DM density profiles, 
find densities that vary by about one order of magnitude [25]. This translates into the same 
order of magnitude uncertainty on the DMz/ rate, due to the direct proportionality between the 
neutrino signal and the local dark matter density through the capture rate. The uncertainty on 
Pdm can also play a role in the setup of capture/annihilation equilibrium in the Earth, giving an 
additional reduction effect. 

An additional astrophysical uncertainty comes from the local DM velocity distribution func- 
tion. Since capture is driven by the relation between the DM velocity and the escape velocity 
of the capturing body (11 km/sec and 620 km/sec at the surface of the Earth and the Sun, re- 
spectively), the high-velocity tail of the DM velocity distribution function may play a role. In 
the case of the Earth, the actual motion of DM particle in the solar system is another relevant 
ingredient which can alter significantly the predicted capture rate and therefore the predicted 
DMz/ rate. Recently this issue has been re-evaluated in [23] , where it has been shown that in the 
Earth the capture rate of DM particles heavier than a few hundreds of GeV may be considerably 
reduced. 

On the contrary, neutrino spectra can be considered as virtually free from astrophysical un- 
certainties. The shape of the spectra depends on the type of particle produced in the annihilation 
process and on its subsequent energy-loss processes (remember that annihilation occurs in a 
medium, not in vacuum) before decaying into neutrinos. As a consequence of the thermalization 
of the captured DM, the density distribution of DM particles within the Sun or the Earth is 
predicted to be 

R 

n{r) = noexp(-rVi?DM) ^dm = (4) 

where r is the radial coordinate, /3 = 2txG^PqR^ /'iTo, po and Tq are the central density and 
temperature of the body (Sun or Earth) and R is its radius. These astrophysical parameters are 
relatively well known, much better than the above mentioned galactic ones. Numerically: 

f 1.76/ GeV for the Earth, 
^ ~ \ 98.3/ GeV for the Sun. 



This means that the size of the production region of DM neutrinos is ~ SOOkmA/lOO GeV/mDM 
in the Earth and ~ 0.01 i?0^1OO GeV/mDM in the Sun. 

The finite size can affect the spectra in two ways: 1) Different DMz/ originate in regions with 
different densities, so that hadrons may loose different amounts of energy before decaying into 
neutrinos. This, however, is not an important effect because the size of the production region 
is small enough that the matter density can be safely considered as constant where neutrinos 
are produced; 2) Neutrino propagation: while in the Earth the production region has a size 
much smaller than the atmospheric oscillation length, in the Sun the size is instead comparable. 
The resulting coherence between different flavors gets however washed-out by the much longer 
eventual propagation up to the Earth. 

In conclusion, the production regions are small enough that performing the spatial average 
according to eq. (jl]) gives a final total spectrum negligibly different than the one obtained by just 
assuming that all DMz/ are produced at the center of the Earth or of the Sun. We will prove this 
statement in section 13. 2[ after discussing our treatment of neutrino propagation. 
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DMz/ spectra and fundamental parameters Since DM particles inside the Earth or the Sun 
are highly non-relativistic, their annihilations occur almost at rest and the main phenomenological 
parameters that determine DMu spectra are the DM mass rriDM and the BR of the basic channels 
into which DM particles may annihilate, as shown in eq. ([1]): qq, ii, uu, W^W~, ZZ and higgs 
particles or mixed higgs/gauge boson final states O |27]. Besides the direct uu annihilation 
channel, neutrinos originate from the decays of the particles produced in the annihilation. In 
the case of quarks, hadronization will produce hadrons whose subsequent decay may produce 
neutrinos. Also charged leptons, apart from electrons, will produce neutrinos. In the case of 
gauge bosons or higgs particles, their decay will produce again leptons or quarks, which then 
follow the same evolutions just described. 

The basic "building blocks" we need in order to calculate DMz/ fiuxes are therefore the spectra 
produced by the hadronization of quarks and by the decay of charged leptons in the Sun and Earth 
cores. Among leptons, only the r is relevant, since muons are stopped inside the Earth and the 
Sun before they can decay [27], and therefore produce neutrinos of energy below experimental 
thresholds for the signal topologies we will discuss later on (up-going muons, contained muons 
and showers in large area neutrino telescopes). For the present discussion we consider neutrino 
energies above 0.5 GeV. In all the other situations, which involve gauge and higgs bosons, we 
can make use of the basic spectra discussed above and calculate the neutrino spectra by just 
composing properly boosted spectra originated from quarks or r, following the decay chain of 
the relevant annihilation final state particle. The method is briefly sketched in Appendix |X] for 
completeness. 

In this paper we are interested in the discussion of the effect induced by oscillations on the 
neutrino signal. We therefore need to calculate the spectra for all three neutrino flavors. We 
model the hadronization and decay processes by means of a PYTHIA Monte Carlo simulation [28] , 
suitably modifled in order to take into account the relevant energy losses. The neutrino spectra 
which we obtain are presented in numerical form, but we also provide an interpolating function 
for all the quark flavors and the r lepton. 

We will not consider effects on the neutrino fluxes arising from the spin of the DM particle. 
In general, the DM spin may control the polarization of primary particles produced in the an- 
nihilation. For instance, if the DM is a Majorana fermion (such as the neutralino) it can only 
decay into tlTr + h.c. (with amplitude proportional to the r mass) while a scalar can decay into 
t^tl and Tjifji with different branching ratios. Only if the branching ratios are the same the 
DMu spectrum is equal to the Majorana case. We will assume that this is the case, studying a 
single rf channel rather than two slightly different tlTl and TRffi channels. Furthermore, when 
discussing direct annihilation into neutrinos (possible for a scalar DM) we will assume that the 
flux is equally divided among the three flavors. 

We now discuss the calculation of the spectra for the relevant flnal states and their distinctive 
features. In this paper we will not focus on a speciflc DM candidate, rather we will attempt a 
more general phenomenological analysis. Our results can therefore be used for any DM candidate, 
by using the basic spectra of primary annihilation particles given here. The full spectrum for a 
speciflc candidate in a speciflc model is then easily constructed by summing up these building 
blocks implemented by the information on the annihilation branching ratios BR/ in that model. 
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2.2 Annihilation into light fermions 

The direct DM DM — )■ vv channel (if allowed with a reasonable branching ratio) usually gives the 
dominant contribution to DMz/ signals: its spectrum is a line at Ey = moM so it gives the neutrinos 
with highest multiplicity and energy. If the DM is a Majorana fermion the DM DM — )■ // 
annihilation amplitude is proportional to rrif, so that the channel is irrelevant and the most 
important fermions are the heaviest ones: bb, rf cc and, if kinematically accessible, ti {i.e. if 
> ^t)- Even in the context of SUSY models the relative weight of their branching ratios 
should be considered as a free parameter: significant deviations from the qualitative expectation 
cr(DM DM — )■ 6fe)/(T(DM DM — )■ rf) = 3m^/m^, which exactly holds in the case of a dominant 
higgs exchange, can arise if staus are much lighter than sbottoms. 

Once a quark is produced, it will hadronize and produce a large number of mesons and 
baryons, which will then decay and eventually produce neutrinos. We calculate the Vei '^n and 
Ur fluxes originated by quark hadronization and lepton decay in the medium by following and 
properly adapting the method of [2^. We improve on previous analyses [5] by calculating full 
spectra for all neutrino flavors (usually only i/^ were considered since the main signal is upgoing 
muons and oscillations have been neglected, except in a few seminal cases [IS])- We also provide 
here neutrino fluxes coming from light quarks: their contribution is usually neglected since they 
mostly hadronize into pions, which are stopped in the medium and do not produce neutrinos in 
an interesting energy range. We show below that for relatively large DM masses neutrinos from 
light quarks should be taken into account in a precise computation. Their main contribution 
occurs through the excitation of c quarks in the hadronization process and subsequent decay of c 
mesons. 

For completeness, we include also the case of DM annihilation into gluons, which may be rele- 
vant for some DM candidate. For instance, in the case of neutralinos, gluons can be produced at 
one loop level: even though this channel is usually subdominant, it can provide some contribution 
in specific portions of the SUSY parameter space, especially for light neutralinos. 

Annihilation inside the Earth As previously discussed, the annihilation process occurs pri- 
marily in the center of the Earth, where the density is p = 13 g cm~'^. Therefore the particles 
produced in DM annihilations may undergo energy loss before decay. 

In the case of charged leptons, the energy loss process is calculated by means of the Bethe- 
Bloch equation. The typical stopping time is of the order of Tstop = 2 ■ 10"^° sec. This has to be 
compared to the boosted lifetime 7rdec, where r^ec = 2.2 ■ 10~^ sec for the muon and t^cc = 3 ■ 10~^^ 
sec for the tau. We see that for leptons with energies up to 1 TeV, muons are always stopped 
before their decays, while taus may decay as if they were in vacuum . In order to take into 
account the small deviations from the limit described above we adapted the PYTHIA code to 
allow free lepton decay if 7rdec < Tstop, otherwise the lepton is stopped and then it decays. 

As for the hadrons, the situation is different if the jets are produced by light or heavy quarks. 
The interaction time in a material with density p is [27] 

Tint = [n aintv]~^ = 5 ■ aint/3]"^sec (6) 

where (Jint denotes the typical interaction cross section for a hadron. 
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For hadrons made of light quarks (Tint ~ 20 mbarn [29], which imphes Tint — 2- 10"^*^ sec. Since 
the typical lifetime of vr and K mesons is of the order of 10"*^ s, light hadrons are usually stopped 
before decay, unless they are very relativistic. We implemented this process in the PYTHIA 
code by letting the hadron decay freely when 7rdcc < Tint, otherwise it is stopped. With this 
modification of the code we take into account the actual lifetime of any hadron and the actual 
energy it has in the fragmentation process. When a very energetic light hadron is produced, we 
therefore do not neglect its decay. This situation however is not very frequent and it can occur 
only for very energetic injected jets. 

In the case of heavy hadrons one has (Tint ~ 14 mbarn for a c or 6 meson and (Tint ~ 24 mbarn 
for a c or 6 hadron [29], giving Tint ~ (2 4- 3) ■ 10~^° sec. The typical lifetime for these hadrons 
is Tdec ~ 10~^^ sec, or less. We therefore may assume that they decay before loosing a significant 
part of their energy. We again implemented a modification of the PYTHIA code which is similar 
to the case of leptons and which takes into account the relevant time scales. 

In addition to energy losses, we should also take into account that interaction of hadrons with 
the medium could lead to the production of additional hadrons. For instance, a heavy-hadron 
collision with the medium may produce additional light hadrons. However these additional light 
hadrons of lower energies are easily stopped, as discussed before, and therefore give a negligible 
contribution to the neutrino flux in our relevant energy range from this process. We therefore 
ignore here this possibility, a consistent assumption under our approximations. 

Our results on the neutrino spectra from annihilations in the Earth are shown in fig. [2] as 
dotted lines. Each spectrum refers to the flux of neutrinos for a given qq or rf pair and for 
different values of ttidm, equal to the energy of the primary jet or r. Antineutrinos are not 
summed up and their fluxes are the same as those of neutrinos, since the initial state is neutral 
respect to all quantum numbers. The plots are shown as a function oi x = E/m-£,M, which is 
defined in the interval [0, 1]. The curves start from the x corresponding to the minimal neutrino 
energy that we consider, Ei, = 0.5 GeV. 

The spectra at production of z/g and u^^ are equal, since light hadrons and muons do not 
contribute to the DMz/ fluxes, and since r, c, b produce an equal amount of z/g and z/^. This 
equality would not hold for neutrinos produced from vr or /i. 

We also see that light quarks contribute to the neutrino fluxes, even though light hadrons 
are stopped. This is due to the fact that a. u, d or s quark has a non vanishing probability of 
splitting into a c quark in the fragmentation process and this process is favored for larger energies 
(for details, see [2S] and references therein). We therefore have c hadrons in the outgoing jets 
also when we inject a light quark. The decay of these hadrons produces neutrino fluxes in the 
interesting energy range. We see that at low neutrino energies around 1 GeV and for mDM larger 
than about 500 GeV the contribution coming from light quarks can even be the dominant one. 
This effect was neglected in previous analyses. 

We provide analytical fitted formula for the spectra. We fitted the MC results with the 
following expression, which proved to be suitable: 

9i^) = ~r~ — '^0(1 + o-iUJ + a2w'^ + a^w^ + a^w^ + a5W^)(l — x)^ + Cox'^^(l — xY^ (7) 
ax 

where w = log^oa^- The values of the parameters are shown in tabled] and table [2] for a sample of 
center-of-mass energies mDM of the primary quark or r, and are also available at [30]. The fitted 
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Figure 2: Neutrino spectra at production. Upper half: the fluxes of electron and muon neutrinos, 
for the seven main annihilation channels and for different masses of the parent DM particle 
(different colors). The solid lines apply to the case of the Sun, the dotted of the Earth. In all 
cases, the spectra of antineutrinos are the same as those of neutrinos. Lower half: the same for 
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6.27/6.40 


0/0 


0/0 


0/0 


30 


0.233/0.674 


-7.56/-3.07 


-8.62/-5.38 


-3.09/-3.29 


0.353/-0.691 


0.297/0.0012 


6.35/7.86 


0/0 


0/0 


0/0 


50 


0.513/0.828 


-2.06/-2.93 


-2.90/-5.87 


-1.62/-4.35 


-0.323/-1.38 


-0.0032/-0.154 


6.81/8.84 


0/0 


0/0 


0/0 


70 


0.328/0.962 


-2.66/-3.78 


-0.903/-8.23 


1.82/-6.68 


1.66/-2.40 


0.388/-0.324 


6.55/9.91 


0/0 


0/0 


0/0 


100 


0.685/1.60 


-1.92/-1.70 


-4.19/-4.59 


-3.65/-3.87 


-1.42/-1.39 


-0.208/-0.185 


7.57/11.1 


0/0.160 


0/6.11 


0/12.8 


200 


0.532/1.41 


-2.17/-2.30 


-3.51/-4.73 


-2.49/-3.42 


-0.727/-1.06 


-0.074/-0.120 


7.58/12.9 


0/0.158 


0/6.11 


0/12.8 


300 


0.279/1.72 


-3.52/-2.68 


-2.93/-4.98 


-1.03/-3.28 


0.118/-0.919 


0.072/-0.094 


6.98/16.8 


0/0.153 


0/6.11 


0/12.8 


500 


0.363/1.35 


-3.02/-3.24 


-3.31/-4.89 


-1.61/-2.74 


-0.136/-0.613 


0.029/-0.045 


7.56/17.5 


0/0.153 


0/6.11 


0/12.8 


700 


0.476/2.90 


-1.60/-3.10 


-1.64/-5.64 


-0.663/-3.50 


0.088/-0.909 


0.044/-0.085 


7.77/23.4 


0/0.152 


0/6.11 


0/12.8 


1000 


0.420/2.19 


-2.08/-3.80 


-1.65/-5.49 


-0.195/-2.84 


0.408/-0.584 


0.095/-0.040 


8.01/26.4 


0/0.151 


0/6.11 


0/12.8 


DM annihilations into qq 


10 


0.0024/0.919 


-2.17/1.04 


-4.97/1.73 


-6.24/5.85 


-4.68/6.21 


-1.38/2.00 


16.0/42.2 


0/0.025 


0/2.59 


0/11.2 


30 


0.038/0.871 


-0.632/3.50 


-5.46/4.96 


-7.21/3.52 


-3.64/1.26 


-0.649/0.183 


17.3/20.4 


0/0 


0/0 


0/0 


50 


0.020/0.405 


2.78/3.44 


5.17/4.79 


5.22/3.18 


3.26/1.03 


0.771/0.132 


12.0/20.4 


0/0 


0/0 


0/0 


70 


0.017/0.0057 


2.28/-2.11 


5.65/-0.717 


5.80/-1.42 


3.44/0.398 


0.733/0.208 


14.2/15.5 


0/0 


0/0 


0/0 


100 


0.019/0.017 


2.25/-0.973 


6.88/-2.07 


8.24/-0.842 


5.16/1.31 


1.10/0.494 


14.6/15.3 


0/0 


0/0 


0/0 


200 


0.017/0.012 


2.28/-0.863 


6.86/-1.64 


7.97/-1.60 


5.49/1.65 


1.14/0.561 


14.3/15.4 


0/0 


0/0 


0/0 


300 


0.015/0.041 


2.05/-1.38 


6.28/-2.16 


7.67/0.967 


6.01/2.60 


1.23/0.635 


14.3/17.3 


0/0 


0/0 


0/0 


500 


0.065/0.047 


3.44/-1.73 


7.82/-1.49 


8.11/2.53 


4.25/3.41 


0.708/0.747 


15.5/19.7 


0/0 


0/0 


0/0 


700 


0.093/0.043 


3.57/-1.50 


7.95/-1.57 


8.14/2.47 


4.04/3.63 


0.638/0.754 


15.7/19.7 


0/0 


0/0 


0/0 


1000 


0.061/0.050 


2.90/-2.07 


8.48/-2.74 


10.2/1.18 


5.57/2.90 


0.873/0.579 


16.3/23.5 


0/0 


0/0 


0/0 


DM annihilations into gluons 


10 


0.050/0.017 


-0.286/0.645 


-1.43/7.23 


-2.12/1.20 


-0.727/-4.95 


-0.011/-2.21 


9.31/9.15 


0/0 


0/0 


0/0 


30 


0.042/0.012 


-2.49/0.394 


-0.522/8.54 


3.44/-2.01 


5.60/-1.81 


2.00/0.043 


8.88/7.62 


0/0 


0/0 


0/0 


50 


0.011/0.802 


-4.03/0.211 


3.80/-2.76 


-0.409/-3.88 


4.15/-1.88 


1.90/-0.310 


7.10/14.5 


0/0 


0/0 


0/0 


70 


0.013/0.532 


-0.695/0.030 


3.66/-2.72 


-3.72/-3.74 


2.55/-1.73 


1.56/-0.271 


6.06/13.8 


0/0 


0/0 


0/0 


100 


0.353/1.01 


-2.00/-0.413 


-6.65/-3.37 


-6.74/-3.66 


-2.61/-1.47 


-0.352/-0.205 


13.1/17.4 


0/0 


0/0 


0/0 


200 


0.082/0.555 


-4.58/-0.850 


-6.82/-3.38 


-5.00/-3.39 


-0.326/-1.18 


0.207/-0.139 


11.1/17.4 


0/0 


0/0 


0/0 


300 


0.052/0.227 


-4.17/-2.12 


-6.56/-3.38 


-5.60/-2.52 


0.399/-0.196 


0.433/0.073 


8.88/17.9 


0/0.090 


0/2.57 


0/8.13 


500 


0.063/0.814 


-5.59/-1.86 


-6.83/-4.34 


-3.22/-3.20 


1.88/-0.802 


0.619/-0.062 


10.4/23.9 


0/1.43 


0/2.41 


0/13.6 


700 


0.069/0.453 


-3.50/-2.42 


0.232/-4.26 


4.58/-2.75 


5.20/-0.358 


1.05/0.023 


11.0/24.0 


0/59.7 


0/4.77 


0/16.4 


1000 


0.235/0.328 


3.74/-2.99 


8.11/-4.19 


7.39/-2.08 


3.63/0.246 


0.561/0.123 


9.14/24.0 


0/481. 


0/5.52 


0/20.6 



Table 1: Fit parameters for the expression g{x) in eq. Give the I'e = i^fi = Ve = spectra at 
production from annihilations in the Earth/Sun. DM masses are in GeV. These data are available 
at J3^. 
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aa 


ai 


02 




a4 


a 5 


b 


CO 


ci 


C2 










DM 


annihilations into 


bb 










10 


5.30/12.9 


-7.48/-4.22 


-29.0/-2a.8 


-37.9/-30.1 


-21.7/-18.8 


-4.69/-4.36 


12.3/13.6 


0/0 


0/0 


0/0 


30 


0.629/0.927 


-14.8/-10.5 


-35.4/-26.1 


-32.9/-24.7 


-13.8/-10.5 


-2.20/-1.70 


10.2/10.9 


0/0 


0/0 


0/0 


50 


0.387/0.766 


-17.3/-11.5 


-36.5/-26.7 


-30.9/-23.8 


-11.9/-9.59 


-1.75/-1.45 


9.86/11.2 


0/0 


0/0 


0/0 


70 


0.290/0.326 


-18.9/-16.4 


-37.0/-28.3 


-29.7/-19.6 


-10.9/-6.00 


-1.51/-0.668 


9.65/10.7 


0/0 


0/0 


0/0 


100 


0.184/0.284 


-21.6/-16.8 


-37.3/-28.1 


-27.6/-19.5 


-9.41/-6.13 


-1.21/-0.715 


9.08/10.8 


0/0 


0/0 


0/0 


200 


0.143/0.381 


-23.9/-15.8 


-37.4/-28.1 


-25.8/-2a.3 


-8.10/-6.60 


-0.959/-0.806 


9.05/12.4 


0/0 


0/0 


0/0 


300 


0.116/0.103 


-26.0/-26.0 


-37.4/-25.7 


-24.0/-11.6 


-6.99/-2.21 


-0.760/-0.129 


9.01/12.3 


0/0 


0/0 


0/0 


500 


0.090/0.190 


-28.6/-20.9 


-36.8/-27.6 


-21.7/-16.3 


-5.71/-4.33 


-0.551/-0.430 


9.00/13.9 


0/980. 


0/5.78 


0/15.3 


700 


0.073/0.085 


-30.7/-28.9 


-36.1/-23.6 


-20.0/-8.85 


-4.84/-1.06 


-0.419/0.030 


8.84/14.1 


0/45.2910 


0/6.72 


0/22.0 


1000 


0.087/0.086 


-29.3/-35.8 


-36.9/-16.3 


-21.2/— 0.119 


-5.32/1.85 


-0.483/0.344 


9.37/23.1 


0/9.00 


0/1.78 


0/15.3 










DM 


annihilations into 


rr 










10 


3.75/3.90 


2.45/2.58 


5.73/6.14 


6.85/7.49 


4.01/4.49 


0.915/1.05 


1.20/1.22 


0/0 


0/0 


0/0 


30 


2.34/2.57 


0.855/1.20 


1.79/2.52 


1.82/2.60 


0.867/1.26 


0.156/0.233 


0.996/1.04 


0/0 


0/0 


0/0 


50 


4.49/4.28 


2.09/2.00 


3.09/2.93 


2.34/2.22 


0.892/0.848 


0.134/0.128 


1.28/1.25 


0/0 


0/0 


0/0 


70 


4.13/4.05 


1.87/1.81 


2.59/2.45 


1.84/1.70 


0.654/0.593 


0.092/0.082 


1.24/1.23 


0/0 


0/0 


0/0 


100 


3.99/4.02 


1.78/1.80 


2.38/2.42 


1.62/1.66 


0.552/0.568 


0.074/0.077 


1.22/1.23 


0/0 


0/0 


0/0 


200 


3.44/3.43 


1.40/1.39 


1.62/1.62 


0.943/0.943 


0.272/0.274 


0.031/0.031 


1.15/1.15 


0/0 


0/0 


0/0 


300 


3.24/3.21 


1.25/1.23 


1.38/1.34 


0.749/0.723 


0.202/0.193 


0.021/0.020 


1.12/1.12 


0/0 


0/0 


0/0 


500 


3.20/3.25 


1.22/1.25 


1.34/1.38 


0.722/0.751 


0.193/0.202 


0.020/0.021 


1.12/1.12 


0/0 


0/0 


0/0 


700 


2.96/3.50 


1.05/1.42 


1.07/1.65 


0.540/0.944 


0.136/0.265 


0.013/0.029 


1.08/1.17 


0/0 


0/0 


0/0 


1000 


2.97/3.00 


1.06/1.08 


1.13/1.14 


0.586/0.588 


0.151/0.150 


0.015/0.015 


1.09/1.09 


0/0 


0/0 


0/0 










DM 


annihilations into 


cc 










10 


-9.25/0.434 


1.33/0.371 


-3.55/-8.04 


-8.46/-16.3 


-6.08/-11.7 


-1.48/-2.93 


32.0/10.3 


0.014/33.7610 


0.348/12.7 


4.75/11.4 


30 


-0.681/0.064 


1.79/-0.626 


-0.550/-9.81 


-2.61/-16.2 


-1.66/-9.84 


-0.328/-2.05 


21.9/7.18 


0.010/0 


0.200/0 


4.90/0 


50 


-2.19/0.012 


3.30/-0.739 


4.15/-9.93 


2.55/-16.5 


0.767/-8.97 


0.091/-1.59 


22.2/4.11 


0.011/0 


0.245/0 


5.14/0 


70 


-0.087/0.652 


-0.198/3.83 


-5.46/5.72 


-6.28/4.05 


-2.72/1.39 


-0.414/0.186 


14.1/9.87 


0.051/0 


0.966/0 


6.81/0 


100 


-4.27/1.05 


3.33/3.67 


4.24/5.23 


2.63/3.57 


0.799/1.18 


0.095/0.151 


26.6/11.1 


0.015/0 


0.307/0 


5.71/0 


200 


0.015/0.031 


1.59/0.326 


-2.99/-4.78 


-6.99/-7.54 


-3.50/-3.60 


-0.546/-0.566 


7.48/5.94 


0.0073/0 


-0.067/0 


5.20/0 


300 


-0.027/0.0075 


0.630/0.650 


-2.12/-5.05 


-2.18/-7.97 


-1.10/-1.92 


-0.194/-0.026 


8.78/6.00 


0.022/232. 


0.512/4.74 


6.42/18.9 


500 


0.323/0.022 


4.35/-1.40 


7.00/-8.57 


5.33/-10.1 


1.98/-3.80 


0.297/-0.483 


10.5/9.25 


0.00083/0.069 


-1.10/2.02 


3.13/6.79 


700 


0.0066/0.013 


0.048/-2.95 


-3.54/-9.83 


-3.66/-8.34 


0.714/-1.50 


0.352/0.017 


7.27/8.55 


0.012/0.561 


0.169/2.66 


6.12/10.4 


1000 


0.0083/0.035 


-2.01/-2.47 


-1.13/-9.60 


2.26/-8.66 


3.59/-2.50 


0.748/-0.240 


13.5/8.54 


0.010/0 


0.083/0 


6.05/0 










DM 


annihilations into 












10 


0.00084/0.012 


1.76/2.20 


0.493/2.52 


-0.811/1.60 


-0.585/0.527 


-0.108/0.066 


18.8/21.2 


0/0 


0/0 


0/0 


30 


0.0032/0.012 


1.06/2.20 


-0.294/2.52 


-1.33/1.60 


-0.723/0.527 


-0.114/0.066 


16.4/21.2 


0/0 


0/0 


0/0 


50 


0.0049/0.042 


2.79/2.45 


3.19/2.57 


1.06/1.38 


0.042/0.429 


-0.017/0.066 


14.8/21.2 


0/0 


0/0 


0/0 


70 


0.0053/0.071 


2.49/2.46 


3.28/2.56 


1.57/1.38 


0.469/0.431 


0.081/0.064 


14.8/21.2 


0/0 


0/0 


0/0 


100 


0.0094/0.217 


1.42/2.72 


-0.133/2.71 


-1.71/1.12 


-0.834/0.156 


-0.111/-0.0047 


14.9/21.3 


0/0 


0/0 


0/0 


200 


0.0066/0.036 


1.33/2.66 


-0.235/2.66 


-1.74/1.14 


-0.306/0.362 


0.048/0.063 


14.9/21.3 


0/0 


0/0 


0/0 


300 


0.044/0.088 


0.368/-3.12 


-1.68/-8.93 


-2.04/-7.62 


-0.696/-2.63 


-0.076/-0.328 


23.0/32.8 


0/0 


0/0 


0/0 


500 


0.012/0.026 


-0.660/-8.28 


-1.80/-8.91 


-1.04/-3.36 


0.566/-0.0037 


0.201/0.128 


23.0/66.9 


0/1.50 


0/1.39 


0/34.0 


700 


0.014/0.080 


-0.566/-7.99 


-1.84/-9.32 


-1.09/-2.99 


0.570/0.129 


0.186/0.121 


23.0/66.9 


0/1.46 


0/1.62 


0/34.0 


1000 


0.114/0.018 


0.00048/-8.85 


-1.94/-7.69 


-1.58/-4.58 


-0.342/-0.993 


-0.018/-0.080 


23.1/67.0 


0/1.48 


0/1.75 


0/34.0 










DM annihilations into g 


uons 










10 


0.054/0.057 


-1.40/-1.05 


-5.22/-5.92 


-4.66/-7.03 


-1.73/-3.33 


-0.235/-0.498 


19.0/18.7 


0/4.65 


0/6.26 


0/12.8 


30 


0.131/0.607 


0.249/-0.739 


-2.65/-5.88 


-4.08/-7.25 


-2.17/-3.56 


-0.392/-0.630 


15.4/21.8 


0/0.719 


0/3.64 


0/12.7 


50 


0.236/0.590 


-0.0054/-1.37 


-3.05/-6.41 


-3.97/-6.82 


-1.89/-2.97 


-0.312/-0.467 


17.1/23.5 


0/0.394 


0/2.27 


0/12.7 


70 


0.568/0.638 


1.48/-1.47 


0.407/-6.39 


-0.625/-6.55 


-0.440/-2.74 


-0.081/-0.416 


18.4/23.9 


0/0.029 


0/0.431 


0/12.7 


100 


0.696/0.568 


1.74/-1.28 


1.03/-5.89 


0.0057/-6.04 


-0.154/-2.49 


-0.033/-0.367 


18.7/22.7 


0/0.245 


0/2.25 


0/12.7 


200 


0.750/0.118 


2.03/-3.84 


1.75/-6.39 


0.682/-3.73 


0.139/-0.686 


0.014/-0.0021 


17.8/22.3 


0/0.198 


0/1.70 


0/12.7 


300 


0.094/0.420 


-3.76/-2.79 


-8.88/-6.76 


-7.42/-5.10 


-2.38/-1.55 


-0.262/-0.168 


16.1/27.5 


0/0.110 


0/1.28 


0/12.7 


500 


0.480/0.160 


0.831/-4.43 


-0.339/-6.73 


-0.728/-3.55 


-0.240/-0.547 


-0.023/0.0066 


18.9/27.5 


0/0.201 


0/1.52 


0/12.7 


700 


0.123/0.112 


-4.50/-4.53 


-9.65/-6.57 


-7.07/-3.70 


-1.90/-0.582 


-0.172/0.0030 


16.8/27.5 


0/0.058 


0/0.698 


0/12.7 


1000 


0.596/0.071 


1.80/6.05 


1.48/12.5 


0.613/8.52 


0.198/3.03 


0.028/0.387 


16.9/24.3 


0/0.081 


0/0.750 


0/12.7 




Table 2: 


Fit paramet 


ers for the 


expression g{x) in eq. 


(0). Give the 


Vr = i^T 


spectra at production 



from annihilations in the Earth/Sun. DM masses are in GeV. These data are available at [30^ . 
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functions reproduce the MC result at a level better than a few percent in all the relevant energy 
range, from 0.5 GeV up to ttidm- The functions g{x) should not be used outside this range. 

Annihilation inside the Sun. The density of the core of the Sun is p = 140 g cm~^, about 
10 times larger than in the Earth, so that energy loss processes are more important than in the 
Earth case. 

In the case of charged leptons, the stopping time is now Tstop ~ 10^^^ sec. Our modification 
of the PYTHIA code takes into account this situation, as described previously. 

The situation for the light-quark hadrons is similar to the case of the Earth: they are mostly 
stopped and therefore they do not produce neutrinos in the energy range of interest. In the case 
of hadrons made by heavy quarks, the situation is now more subtle [27] . Their typical interaction 
time gets reduced by an order of magnitude: Tint ~ (2 3) ■ 10~^^ sec, and becomes comparable 
to the typical heavy-hadron lifetime Tdcc ~ 10~^^ sec (some hadrons decay faster). We must now 
be careful, since these hadrons may loose a fraction of their energy before decaying. In order to 
take into account this effect, we follow [27] where the average energy loss of a heavy hadron in a 
dense medium was studied. For details about the analysis, we refer to [27]. Here we just recall 
the relevant results, implemented in our analysis. 

A c or 6 hadron of initial energy Eo after energy losses emerges with an average energy: 

(E) = E, exp{xo) E,{xo) (8) 
where Ec = Mhadron Tstop/Tdcc, a^o = Ec/Eq and the function Ei is defined as: 

poo g-x 

E.ixo) = / dx. (9) 

The quantity Tstop is defined as Tstop = Tint/ i)- — Z) where Z = Xi Zj and Tj^t is the interaction time 
defined in eq. (E]). The quantity Xi denotes the ratio between the quark and the hadron mass: 
Xi = mj/Mhadron, for i = c,b and Zj = 0.6 for a c hadron and zj = 0.8 for a b hadron. 

We modified the PYTHIA code in order to take into account the energy loss discussed above: 
when a c or b hadron is produced, we first reduce its energy according to eq. (|8]), and then it is 
propagated and decayed by the PYTHIA routines. 

Our results for the neutrino spectra from annihilations in the Sun are shown in fig. [2] as solid 
lines. We see that the spectra are a little softer than in the Earth case, due to hadron energy 
losses. The effect is more pronounced for larger center-of-mass energies, since in this case hadrons 
loose a larger fraction of their initial energy. Also in this case we fitted the distributions with 
the same fitting formula of eq. ([7]), and reported the parameters in tabled] and table |2J They are 
again available at [30] . 

2.3 Annihilations into W+W and ZZ 

The lifetime of W, Z gauge bosons is short enough that their energy losses can be neglected. They 
therefore decay into quarks and leptons as in vacuum, but then their decay products hadronize 
and decay, loosing energies as discussed in the previous paragraphs. We can therefore calculate 
the neutrino spectra by applying the results for quarks and leptons and by using the formulae 
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given in Appendix |Al To the resulting spectra we than have to add the production of 'prompt' 
neutrinos by the decays W — z/^f and Z — uiUi, that give neutrino hues in the reference frame of 
the decaying boson. When the boson is produced with an energy Eb, the neutrino hne is boosted 
to a flat spectrum: 

dN BRi , 1-/3 1 + /3 

— — = — — m the range — - — < a; < — - — (10) 
ax p 2 2 

where BRj is the branching ratio for the prompt decay of the gauge boson and (5 is velocity of 
the gauge boson. As a check to our calculation, we produced a few sample cases of neutrino 
spectra from W and Z with the PYTHIA code and compared them to our analytical results. The 
agreement is well under the MC statistical error. 

Our results are shown in fig. |2]as dotted lines for the Earth and as solid lines for the Sun. DM 
annihilations into vector bosons produce a harder DMz/ spectrum as compared to annihilations 
into rf, 66, cc, gg, thanks to prompt neutrino production. This is a dominant feature in the 
spectrum as long as Eb = ?Ti£)M is not too much larger than M^y, Mz- 



2.4 Annihilation into tt 

The lifetime of the top quark is extremely short too, which allows us to consider it decaying before 
any energy loss is operative. Also in this case we build the spectra for the tt case as described 
in Appendix |Al by using the decay chain: t h + W followed by W decay, as discussed in the 
previous paragraph. 

Notice that we consider a pure SM decay for the top quark. In two-higgs doublet models 
like e.g in supersymmetric extensions of the SM, there may be additional final states for the top 
decay, due to the presence of a charged higgs: t — )■ 6 + , followed by h hadronization and 
decay. Similarly, we do not consider DM annihilations into new particles, like e.g. H^H^ . 



2.5 Annihilation into higgs bosons or higgs+gauge bosons 

DMz/ can also be generated by channels involving higgs particles in the annihilation final state. We 
can safely ignore energy losses also for the higgses and directly apply the method of Appendix \^ 
We do not explicitely provide results for this case, because even within the SM the Higgs decays 
remain significantly uncertain until the Higgs mass is unknown. Furthermore, Higgs decays can 
be affected by new physics: e.g. in SUSY models the tree-level Higgs/fermions couplings differ 
from their SM values. 

Higgs decays do not produce prompt neutrinos (because of the extremely small neutrino 
masses) so that only soft neutrinos are generated, even softer that in the qq case. Whenever a 
higgs is produced in conjunction with a gauge boson prompt neutrinos from the gauge boson will 
be present. 



3 Neutrino propagation: oscillations, scatterings,... 

We need to follow the contemporary effect on the neutrino fluxes from DM annihilations (presented 
in the previous section) of coherent flavor oscillations and of interactions with matter. 
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The appropriate formalism for this, that marries in a quantum-mechanically consistent way 
these two aspects, consists in studying the spatial evolution of the 3x3 matrix of densities 
of neutrinos, p{E^), and of anti-neutrinos, p{E^). We will indicate matrices in bold-face and 
use the flavor basis. The diagonal entries of the density matrix represent the population of 
the corresponding flavors, whereas the off-diagonal entries quantify the quantum superposition 
of flavors. Matrix densities are necessary because scatterings damp such coherencies, so that 
neutrinos are not in a pure state. The formalism is readapted from [21], where it was developed 
for studying neutrinos in the early universe. 

The evolution equation, to be evolved from the production point to the detector, has the form 



dp _ j^ijj p\_^ 
dr ' dr 



_l_ dp 

cc dr 



dp 

NC dr 



'ID 



with an analogous equation for p. The first term describes oscillations in vacuum or in matter. 
The second and the third term describe the absorption and re-emission due to CC and NC 
scatterings, in particular including the effect of Vr regeneration. The last term represent the 
neutrino injection due to the annihilation of DM particles. The average over the size of the 
production region can be approximately performed as described below. Note that there is no 
neutrino-neutrino effect (i.e. the evolution equation is linear in p) because neutrino fluxes are weak 
enough that they negligibly modify the surrounding environment. In particular Pauli blocking 
(namely: the suppression of neutrino production that occurs due to fermion statistics if the 
environment is already neutrino- dense), important in the early universe and in supernovae, can 
here be neglected. 

We will discuss each term in detail in the following sections. 

In the case of neutrinos from the center of the Earth, the formalism simplifies: indeed, neutrino 
interactions with Earth matter only become relevant above 10 TeV. Since typical DM particles 
have the correct abundance for m-Qu ^ TeV, we can ignore such interactions in the Earth and 
only oscillations need to be followed. Moreover, taking into account that the initial spectra $° 
do not distinguish e from /i, $e = '^'J^ = (as discussed in sec. [2]), and that the oscillation 
probabilities obey Pre + -Pt/^ + Ptt = 1? the oscillated fluxes are given by 

$,(E,) = ^1^{E,) + Pr,{E,) [^1{E,) - ^1^{E,)] I = {e, /i, r} (12) 

An analogous result holds for anti-neutrinos. P^-^, the conversion probability of a Vr into a 
neutrino of flavor is easily computed with the standard oscillation formalism described below 
and is plotted in fig. [T|d. 



3.1 Oscillations 



Oscillations are computed including the vacuum mixing and the MSW matter effect 
effective Hamiltonian reads 

TY) ' Tti > r 

H = + v^Gf [Ne diag (1, 0, 0) - ^ diag (1, 1, 1) 



The 



(13) 
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where m is the 3x3 neutrino mass matrix. One has m)m = V ■ diag(m^, mg, m|) ■ where 
^^1,2,3 > are the neutrino masses and V is the neutrino mixing matrix. We define the solar 
mixing angle as tan^sun = IK1/K2I5 the atmospheric mixing angle as tan^atm = l^^s/K-sl and 
sin ^13 = |V"e3|. Ne{r) and Nn{r) are the number density of electrons and neutrons in the matter, 
as predicted by solar and Earth models [331 El] • The above Hamiltonian applies to neutrinos; for 
anti-neutrinos one has to replace m with its complex conjugate and fiip the sign of the MSW 
term. The difference between the matter potential for and Uj. that arises only at one loop 
order, becomes relevant only at E^, > fewTeV so we can neglect it. Finally, notice that matter 
effects suppress oscillations of Ue, since they encounter no MSW level crossings. 

In the following we assume the present best fit values for the mixing parameters (from [36] ) 

tan' e^un = 0.45, ^atm = 45°, 

Aml^^ = 8.0 10~^eV^ |Am^^^| = 2.5 10"^ eV'. 

We assume 613 = and we will later comment on how a non-zero ^13 would marginally modify 
our results. 



3.2 Average over the production region 

Neutrinos are produced in the core of the body (Earth or Sun) over a region of size -Rj3M) as 
discussed in section [2l so in principle their propagation baseline is different depending on where 
they originate. However, since Rbm turns out to be smaller than the size of the object, to a 
good approximation one can take into account oscillation effects assuming that all neutrinos are 
produced at the center of the production region, with the following effective density matriji0 

ipfh= E (f^-f^^™)(Po).Ye-«BM(^^™-^^.)V8(f;*^t;^.,^) (14) 

m,n,i',j' 

where po is the diagonal matrix of the total initial fluxes po{E^) = diag ($g, $°). U and H 
are the energy-dependent neutrino mixing matrix and hamiltonian eigenvalue at the center of the 
body, to be computed diagonalizing the Hamiltonian in matter, eq. ( !T3|) . 

For a better intuitive understanding of the physical meaning of eq. f|T^ one can neglect the 
small oscillation effects driven by 613 and Am'^^^ (since the oscillation length of the latter is much 
larger than the production region) and keep only the oscillations driven by Amatm, thus reducing 
the oscillation to a 'V^ o Uj." case. Such oscillations are not affected by matter effects so that 
U = V and the eigenvalues if2,3 reduce to m^^/Ey. Then the effect of the exponential factor in 
eq. f|T^ is to damp the coherence between the z/2 and z/3 mass eigenstates, i.e. the off-diagonal 
elements of the density matrix (which express the superposition of different states) are suppressed. 
In other words, in the limit of a complete damping the effective density matrix is diagonal and 

^We sketch the proof. To leading order in i?DM the distribution of neutrinos as a function of their path- length 
i + (5 is n{5) cx exp(— 2(5^/_Rqj^), where L is the distance of the detection point from the center of the body and 
6 spans the production region. The factor of 2 accounts for the two DM particles in the annihilation initial state. 
Oscillations can be decomposed as C/(L + 5) = U{L) ■ U{S), with U the time evolution operator. Averaging over 6 
gives eq. p^ . 
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Figure 3: Energy distributions of v (red) and v (blue) produced with energy El by one NC DIS 
interaction of a ^V^ with energy E^. The energy is plotted in units of E^. Continuous line: in 
normal matter, where Np ^ Nn- Dotted line: around the center of the Sun, where Np ^ 2A^„. 

composed of an average of the initial z/^ and z/r fluxes. This is exactly the case for neutrinos with 
very small E^. At larger energies, the damping effect is milder and indeed one can follow the fast 
oscillations. In fig. [T|d the result is well visible. In short: the spatial average over the slightly 
different baselines produces some partial flavor equilibration and some loss of coherence. 

3.3 NC scatterings 

NC scatterings ^vW <H- '•uW effectively remove a neutrino from the flux and re-inject it with a 
lower energy. So they contribute to the evolution equation as: 

dEl^{E^,El)p{E^) + dEl^{El,E^)piEl) (15) 

where 

Fnc(^., El) = Np{r) diaga{uep -> u'^X) + iV„(r) diag(T(z/^r2 ^ u'^X) (16) 

The first term describes the absorption: the integral over El just gives the total NC cross section. 
The second term describes the reinjection of lower energy neutrinos: their spectrum is plotted in 
fig. [3l We see that it negligibly depends on the chemical composition Np/Nn, where Np and Nn 
are the number densities of protons and neutrons. We use the iVe(r) = Np{r) and iV„(r) profiles 
predicted by solar and Earth models [33| IM] . In the Sun Np/Nn varies from the BBN value, 
Xp/Nn ~ 7 present in the outer region r/i?0> 0.3, down to Np/Nn ~ 2 in the central region 
composed of burnt ''He. The Earth is mostly composed by heavy nuclei, so that Np and Nn are 
roughly equal. 

3.4 CC absorptions and Ur regeneration 

The effect of CC interactions to the evolution of the neutrino fiuxes can be intuitively pictured as 
follows. The deep inelastic CC process on a nucleon {^vX — ^J^'X) effectively removes a neutrino 
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Figure 4: Energy distributions of neutrinos regenerated by CC scatterings of a ^v'j. of given energy 
Ey^, produced by one ^Vi/nucleon scattering. The blue upper curves are fr-^riE,^^, E'^), and the 
red lower curves are fT^e,^iiE^^, E^), plotted for several values of the incident energy E^^. 



from the flux and produces an almost collinear charged lepton. The produced by decays 
promptly, before loosing a significant part of its energy in interactions with the surrounding 
matter, and therefore re- injects secondary fluxes of energetic neutrinos 



— e~ Z/g Vt 
IT Uf, Vr 



XUr 
— 7- e+ Ue Vt 
Uf, Ur 



;i7) 



with branching ratios BRx = 64.8%, BRg = 17.84% and BR^ = 17.36% for hadronic, electron 
and muonic decay modes respectively. In this way, besides that always re-appears from 
decays, in ^ 35% of cases also or 'l7}j are produced, enlarging the total neutrino flux that 
reaches a detector. Note that e, /i and hadrons produced by CC scattering and by r decays loose 
essentially all their energy in the matter and are absorbed or decay into neutrinos with negligibly 
small energy for our purposes. 

The CC contribution to the evolution equation of the density matrices is therefore 



dp 

dr 



dp 

dr 



{FccP} 



CC 



dEl 



2 / 

+ne_^p^^(£'™)rcc(-£^™)/-f^e,^j(-E'", Ey) 

{fccP} (dEil 



CC 



2 / 

-fne,^Prr(-E^'')rcc(-E™)/r^e,^(-E'", Ey) 



fl8a) 



:i8b) 



The first terms describe the absorption; their anticommutator arises because loss terms correspond 
to an anti-hermitian effective Hamiltonian such that the usual commutator gets replaced by an 
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anticommutator (see the full formalism in [3T]). The second terms describe the 'z/,- regeneration'. 
We explicitely wrote the equations for neutrinos and for anti-neutrinos because they are coupled 
by the second terms. 

In the formulae above, II^ is the projector on the flavor vi. e.g. 11,- = diag (0, 0, 1). The Fee, 
Fee matrices express the rates of absorption due to the CC scatterings and are given by 



Fcc(^.) = diag (rcc, r^c. r^c), r^c = ^vi^) ^{^iV ^ m + iVn(r) a{v,n ^ (19a) 

fcc(i?.) = diag(r^c,f^c'r^c), r^cc = ^^(0 <y{^tv ^ l^) + ^nir) a{m ^ Ix). (i9b) 



Deep inelastic scatterings of 'z/' on nucleons are the dominant process at the energies involved 
{Ey ^ GeV), so the cross sections a{^vW — )■ ^I'X) (reported in Appendix [B]) are the only ones 
relevant. Scatterings on electrons have a cross section which is ~ iTie/mN smaller and would 
become relevant only at energies E^, > Te V [39] . 

Notice that the matrix Fee is not proportional to the unit matrix because at the relevant 
neutrino energies the 'l7^iV cross sections [IQ] are suppressed with respect to the corresponding 
^''^e,fj,N cross sections by the kinematical effect of the r mass. E.g. at E^, = 100 GeV gives 
a 30% suppression. In particular, this implies that the coherence among '•vl and 'l^e,^ is broken 
by the CC interactions and the formalism is taking this into account. A non trivial consequence 
(interactions increase the oscillation length) is discussed in appendix O 

The functions /{E^, E'^) are the energy distributions of secondary neutrinos produced by a CC 
scattering of an initial neutrino with energy E^^ . They have been precisely computed numerically 
as described in [38]. In the computation of r decay spectra we have taken into account the 
sizable widths of final state hadrons, which produce a significant smearing with respect to fig. 10 
of [H] where such widths are neglected. Fig. H] shows our result for the neutrino spectra from Uj. 
regeneration. The integrals of the z/^ curves equal to one, because z/^ are completely regenerated, 
with lower energy (the curves are peaked at small E^/E^^). The integrals of the z/g,^ curves have 
a value smaller than one, equal to the branching ratio of leptonic r decays. The curves depend, 
but quite mildly, on the incident neutrino energy E^^, mainly due to the finite value of m^: 
neutrinos with lower energy loose a smaller fraction of their energy, because the energy stored in 
the r mass becomes more important at lower energy. We assumed that r~ and r"*" have exact 
left and right helicity respectively; this approximation fails at energies comparable to (say 
El,, Er < 20 GeV [10]), where absorption and regeneration due to CC scatterings become anyhow 
negligible. 

The / functions do not significantly depend on the chemical composition: in the plot we 
assumed Np/Nn = 2 as appropriate in the center of the Sun. Writing x = E^/E^^, these functions 
are normalized to the branching ratios of decays given above: 



Given the ingredients above, it should now be apparent how the second terms in eq.s f|T8|) 
incorporate the CC processes. In words, focussing for definiteness on the case of neutrinos (an- 
tineutrinos follow straightforwardly): a z/^ of energy i?™, described by the density prr, interacts 
with a rate Tq^ and produces secondary z/^, z/g and with energy E^, that contribute to the 
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corresponding diagonal entries of the density matrices p and p. Integrating over i?™ gives the 
total regeneration contribution. 

4 Neutrinos from DM annihilations in the Earth 

In this section we show the results concerning the signal from DM annihilations around the center 
of the Earth: the energy spectra at detector of neutrinos and antineutrinos of all flavors and the 
energy spectra of the main classes of events that they produce. 

Fig. [5] displays the neutrino spectra dN/dE^, from the main annihilation channels DM DM 
— 7- bb, T^T~ , cc, ti, W^W~ , ZZ, normalized to a single DM annihilation]^ A linear combination 
of these basic spectra, weighted according to the BRs predicted by the specific DM model of 
choice and rescaled by the appropriate geometric factors, will give the actual neutrino signal at a 

dN _ Tann V^gp^iV^ _ 0-2 T,,,^ V BR "^^^ (20) 

dtdS dE^ AttRI^ ' dE^ sec ■ m2 lO^Vsec ^ 'dE^' ^ ' 

where the sum runs over the annihilation channels with branching ratios BRj, is the Earth 
radius, and Fann is the total number of DM annihilations per unit time. As already discussed, 
this latter quantity is strongly dependent on the particle physics model under consideration and 
also on astrophysics, and can carry a large uncertainty. When we need to assume a value for it, 
e.g. to compare with the background or with the existing limits, we choose 





AOOGeV\ 


lEarth- 1 









(21) 



In the neutralino case, samplings of the MSSM parameter space find a wide range of 10^~^^ 
annihilations per second, that decreases for increasing rriDM- So our assumption is realistically 
optimistic. 

We show plots for three different values of the DM mass, which give qualitatively different 
results and (in the case of the signal from Earth) well represent the general situation: 

1. rriDu = 50 GeV < Mw,z so that only annihilations into leptons and quarks (other than 
the top) are allowed. Varying m^M in this range the unoscillated fluxes rescale trivially; 
oscillated fluxes also rescale but of course keeping their first dip and peak at fixed energy, 
as described below. 

2. Mwz < JTiDM = 100GeV<mt so that annihilations into vector bosons are kinematically 
allowed, with kinetic energy comparable to their mass. As explained in section [2] this gives 
a characteristic threshold feature: direct decays of W, Z give neutrinos in the energy range 
of eq. (fTOl) (producing the peaks in fig.s|5]), and neutrinos with lower energies are produced 
by secondary decay chains (producing the tails). 

^ These fluxes are available at [50] . 
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Figure 5: Neutrino spectra generated by one DM annihilation around the center of Earth. The 
plots show the spectra of the three neutrino flavors (the three rows) and assume different DM 
masses (the three columns). Each plot shows the open annihilation channels DM DM — )■ bb, t^t~ , 
cc, tt, W^W~ , ZZ . The DM DM — vv channel (not shown) would produce a line at = moM- 
The dotted lines show the spectra without oscillations while solid lines are the final results after 
oscillations. The dashed lines in the upper-left panel have been computed with = 0.1 rad for 
illustration (see text); all the other results assume 9is = 0. Neutrino and anti-neutrino spectra 
are roughly equal: we here show (2$i/ + $p)/3, in view of a{vN) ~ 2a{uN). The shaded region 
is the atmospheric background, normalized relative to DMu as assumed in eq. (12T]) . 
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3. niDM = 200 GeV > rrit so that also annihilations into top quarks are allowed. Since the 
subsequent decay t Mu is a 3-body process, it does not give threshold features. Z, W 
bosons are so energetic that their threshold features are minor. No new notable features 
appear going to higher moM- If the DM is a neutralino only annihilations into W~^W~, ZZ, ti 
(and possibly higgses and SUSY particles) are relevant. 



The atmospheric neutrino background. In all our figures, the shaded region is the back- 
ground of atmospheric neutrinos, computed as predicted by FLUKA [12] (at the SuperKamiokande 
site) and taking into account atmospheric oscillations. The unknown DMi^ signal is compared 
with the known magnitude of the background assuming the annihilation rate in eq. fl2T|) . 

Since the signal comes from the center of the Earth, the background of atmospheric neutrinos 
can be suppressed exploiting directionality: in the figures we applied an energy- dependent cut on 
the zenith-angle, keeping only neutrinos (and, later, events) with incoming direction that deviates 
from the vertical direction by less than 



where E is the energy of the detected particle and irtN ^ GeV is the nucleon mass. Such a choice 
can be understood as follows. First, the finite size of the DM annihilation region implies that the 
signal comes from a characteristic angular opening ~ Rjyu/R® ~ ^/raj^jEy, where the last 
relation makes use of eq. ^ and of the fact that ~ moM/few. Furthermore, the kinematical 
angle ~ 0.30 {m]s[ / E^Y'^^ [13] between the incident neutrino and the produced lepton must be 
taken into account and gives a comparable effect. Finally, to these angles the effect of the angular 
resolution of detectors should be added. For Cerenkov neutrino telescopes under construction 
such as ANTARES, ICECUBE and the future km^ detector in the Mediterranean this resolution 
is 59 < 1°. For AMANDA and Super-Kamiokande the mean angular resolution is 66 ~ 2° or 
larger, hence the angular cut may be larger than what we apply. 

In summary, a more realistic dedicated analysis of the angular (and energy) spectrum will 
be certainly necessary to disentangle in the best possible way the signal from the atmospheric 
background, but our approximation in eq. f l22|) is a reasonable cut applicable to many experiments. 
We stress that the atmospheric background in the small cone around the vertical can be accurately 
and reliably estimated by interpolation of the measured rates in the adjacent angular bins where 
no DMi> signal is present. 



The effect of oscillations. In fig. [S] the dotted lines show the spectra without (i.e. before) 
oscillations: these spectra have been already described in section |2l The final spectra (solid lines) 
are in many cases significantly different. This is also illustrated in fig. [T5k for a few selected cases. 
Oscillations driven by Aml^^ and 6'atm are the main effect at work in fig. [51 They convert 
O 'l^Ji at E^ < 100 GeV (at larger energies the oscillation length is larger than the Earth 
radius) and thus are of the most importance when the initial ^I/}^ fluxes are significantly different 
from the 'l?^ fluxes. This happens e.g. in the case of the DM DM—)- rf annihilation channel: a r 
decay produces one Ur, and just about 0.2 j/^ with little energy; oscillations subsequently convert 
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DM mass 


DM annihilation channels in the Earth/Sun 
vv hh TT cc qq tt ZZ W^W 


50 GeV 
100 GeV 
200 GeV 
400 GeV 
1000 GeV 


1/0.75 0.50/0.67 3.9/3.2 0.32/0.59 0.48/0.66 -/- -/- -/- 
1/0.55 0.70/0.63 2.0/2.7 0.49/0.55 0.45/0.63 -/- 1.0/0.75 1.1/0.75 
1/0.30 0.86/0.55 1.3/1.9 0.75/0.50 0.55/0.58 1.0/0.64 1.0/0.45 1.0/0.47 
1/0.1 0.95/0.44 1.1/0.91 0.91/0.42 0.77/0.51 1.0/0.31 1.0/0.19 1.0/0.20 
1/0.02 0.99/0.32 1.0/0.28 0.98/0.34 0.93/0.42 1.0/0.11 1.0/0.05 1.0/0.06 



Table 3: Ratios of through-going muon rates 'with' over 'without' the effects of the neutrino 
propagation, for DM annihilations around the center of the Earth/Sun. E.g. the bottom-right 
entry means that, formDM = 1000 GeV, the rate is unaffected z/DMDM — )■ W^W~ annihilations 
occur in the Earth, and the rate gets reduced to 0.06 of its value if annihilations occur in the Sun. 

'^T '^fi and significantly enhance the rate of /i events. For instance, neglecting oscillations the 
XX ~^ bb annihilation channel for neutralinos x is regarded as a more significant source of v^ 
than the XX ~^ 'tt' channel, because of the relative branching ratio of 3{mb/mr)'^ (the precise 
value depending on stau and sbottom masses). Oscillations partly compensate this factor as 
quantitatively shown in table |3l that summarizes the relative enhancements or reductions due 
to oscillations, on the rate of through-going muon events (see below) for different annihilation 
channels and for different DM masses ttidm- 

Oscillations also distort the energy spectrum of neutrinos, in the ways precisely shown in 
fig. |5l Table IHreports the mean *z7|j energy after oscillations, for different DM masses and different 
annihilation channels. Notice that, when kinematically open, the ZZ, W^W~ channel remains 
harder than rf, that is usually quoted as the source of a hard neutrino spectrum. 

It is also worth noticing that the DM neutrino signal comes from a distance L ~ R^, while the 
background of up-going atmospheric neutrinos from L ^ 2/?®. Indeed oscillations produce dips 
in the background atmospheric t'^'s at energies ^ ^m1^^R(^/2Tx{n — 1/2) where n = 1, 2, 3, . . . 
{Ey 26 GeV for n = \) and in the background of atmospheric v^s at E^, ^ ^m1^^R(Q/2Tin ^ 
13GeV/n. Since uncertainties on the determination of |Am^^.jjj| from atmospheric experiments 
are still significant, all above energies could have to be rescaled by up to ±30%, so that our results 
would be somewhat affected. 

Oscillations driven by Am^^^ and 6*3^1 have a little effect (at variance of what will happen for 
DM annihilations in the Sun). Finally, let us comment on the small effect of a non vanishing ^13 
on the fluxes. If 6*13 = 0, 'z7e are decoupled from the oscillations driven by Am^^^^ and ^atm, so their 
spectra are not affected: the solid lines (oscillated results) are actually superimposed to the dotted 
ones (no oscillations) in flg.O Choosing instead ^13 = 0.1 rad, we plot for illustration in the upper- 
left panel of flg. [5] the resulting spectra (dashed lines). It is evident that the modiflcations are 
small and mainly concentrated at low energies, where the atmospheric background is dominant. 
This behavior is readily understandable in terms of eq. f lT^ and flg. [H that plots the conversion 
probabilities as function of the energy. With the same tools, one sees that for the other flavors or 
for more energetic neutrinos, the effects of ^13 are even smaller or completely negligible, so we go 
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DM mass 


DM annihilation channels in the Earth/Sun 
uiy hi Tf cc qq ti ZZ W+W' 


50 GeV 
100 GeV 
200 GeV 
400 GeV 
1000 GeV 


100/90 11/11 32/33 12/11 5.7/4.8 / / / 
100/80 11/9.1 25/31 11/8.8 3.2/3.8 -/- 33/29 34/31 
100/62 12/7.3 22/26 12/6.9 2.8/2.8 15/12 33/22 34/24 
100/35 11/5.5 22/19 12/5.0 2.2/2.1 15/8.9 33/15 35/17 
100/9.5 10/2.8 24/9.6 9.9/2.9 3.9/1.3 15/5.0 33/7.0 36/8.0 



Table 4: Average percentage energies in units ofmDu of ^V^^ produced by DM annihilations around 
the center of the Earth/Sun, computed for various annihilation channels and for various DM 
masses. E.g. the bottom-right entry means that DM DM — W~^W~ annihilations with m^M = 
1000 GeV produce 'v^^ with average energy equal to 36%-mDM = 360 GeV if occurring in the Earth 
and to 8% ■ mDM ^/ in the Sun. In the Earth the dependence on ttidm is due to oscillations ( more 
important at lower niD^) and to energy losses of primary particles (more important at higher 
"^DMy'- In the Sun oscillations give sizable effects for any mDM? Ojnd absorption is significant for 
moM > 100 GeV. 

back to the assumption ^13 = in the other panels of fig. |5] and in all other results from now on. 

Let us summarize the impact of oscillations, making reference to fig.s |5] and [T5b: the Ve fiux 
is unchanged; the fiux is significantly increased only for the r annihilation channel; the Vt- flux 
increases in the case of annihilation into h or c. We next compute the energy spectra of the main 
topologies of events that contribute to the measured rates in detectors. 

4.1 Through-going muons 

Through-going /i^ are the events dominantly generated by up-going ^V^^ scattering with the water 
or (more importantly) with the rock below the detector and that run across the detector. Their 
rate and spectrum negligibly depends on the composition of the material: for definiteness we 
consider the rock case. We compute their spectra by considering all the muons produced by 
neutrinos in the rock underneath the detector base, following the energy loss process in the rock 
itself |13] and collecting all /z^ that reach (with a degraded energy) the detector base. We ignore 
through-going muons produced by ''V\ scatterings with the matter below the detector that produce 

that decay into /i^, as these give only a small contribution ( < 10%) to the total rate. 

Fig. |6] shows their expected spectrum in km~^yr~^. Note that, since both the scattering cross 
section and muon path-length are roughly proportional to the neutrino energy, assuming the 
annihilation rate of eq. (12T]) we get a total fiux that roughly does not depend on mDM- Also, 
note that due to the strong increase of the //^ fiux with the energy of the neutrino, annihilation 
channels that produce very energetic neutrinos (such as DM DM—)- vv) give a much larger fiux 
than channels that produce soft neutrinos (such as DM DM—;- hh). Therefore in fig. [6] we had to 
rescale these fiuxes by appropriate factors. 

In Table |3] we present the ratio of the rates of through-going muons with and without oscilla- 
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Figure 6: Spectra of through-going ii^ and fi~ (summed) generated by DM annihilations around 
the center of Earth. The plots assume the DM annihilation rate of eq. ( 12T]) . different DM masses 
and show the main annihilation channels. For better illustration, some channels have been rescaled 
by the indicated factor. The shaded region is the atmospheric background. 
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Figure 7: Spectra of fully contained fi^ and fi'~ (summed) generated by DM annihilations around 
the center of Earth. We assume the DM annihilation rate of eq. (121 p and a detector with Mton-year 
exposure. The vv channel gives a fi^ spectrum peaked at ~ "W-dm? due to the monochromatic 
parent spectrum. 
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Figure 8: (Idealized) energy spectra of showers generated by DM annihilations around the center 
of Earth, in a detector with Mton-year exposure and unit detector efficiency, assuming the DM 
annihilation rate of eq. fl^T]) . 24 



tions, showing how oscillations affect such experimental observable. 

Present data constrain the total /i^ ffux to be below • yr [TT], the constraint being 

stronger when fi^ are more energetic. The ffuxes in fig. E] obey such limits, except probably the 
case of 100% annihilation into the hardest channel uu that was not considered in the AMANDA 
analysis [TT] . 

4.2 Fully contained muons 

Fully contained muons mean that are created inside the detector and that remain inside the 
detector, such that it is possible to measure their initial energy. We compute their energy spectra 
convoluting the neutrino fluxes plotted in fig. [5] with the cross section in the detector volume. 
They are shown in fig. [TJ where the results are normalized assuming the DM annihilation rate of 
eq. fl2Tl) . and considering a detector with active mass times live-time equal to a Mton-year. 

The extent to which muons can fit into this category depends of course on the size and geometry 
of the detector, because more energetic muons travel a longer distance, and on the possibility to 
apply containment requirements. For instance, a km'^ detector in ice or water (mass 1000 Mton) 
contains muons up to about 100 GeV. However detectors of such large sensitive mass are being 
built with the focus on discovery: such sizes inevitably impose to sacrifice the granularity of the 
detector, implying higher energy thresholds and poorer energy resolution. Neutrino telescope 
detectors may have an insufficient granularity of photo-tubes to allow a safe containment cut. 
ANTARES attempted a study [H] of the energy reconstruction from the muon range for contained 
events but the efficiency at sub-TeV energies is affected by luminous backgrounds in the sea. In 
ICECuBE-like detectors a good energy reconstruction (of the order of ±30%) is achieved above 
the TeV range, which leaves small room for WIMP fluxes. 

A better energy resolution would help in discriminating the signal from the atmospheric back- 
ground (concentrated at lower energies) and would allow to study the properties of the signal. For 
example, the smaller SK detector achieved a 2% in the energy resolution of charged particles, and 
measured quite precisely the energy for the single ring contained events, with < 10 GeV [15]. 
At higher energies neutrino collisions are dominated by deep inelastic scattering interactions, that 
produce multiple final state particles making more difficult to tag the event and to measure their 
energies. In principle their total energy is more strongly correlated to the incoming neutrino 
energy; we here compute the energy spectrum of muons only. A water Cerenkov Mton detector 
could isolate fully contained events up to (20 ^ 30) GeV (depending on its geometry) achieving a 
similar energy resolution as SK. 

4.3 Showers 

So far we considered the traditional signals generated by CC ^z7}i scatterings. We now explore the 
shower events generated by: 

(1) CC scatterings of ^v\. We assume that the total energy of the shower is E^^. This is 
a simplistic and optimistic assumption: the appropriate definition is detector-dependent. 
One can hope that showers allow to reconstruct an energy which is more closely related to 
the incoming neutrino energy than what happens in the case of the /i^ energy. 
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(2) NC scatterings of 'i^e,/i,T- At given energy NC cross sections are about 3 times lower than 
CC cross section. We assume that the energy of the shower is equal to the energy of the 
scattered hadrons. 

(3) CC scatterings of 'z7^, generating that promptly decay into hadrons. We assume that the 
shower energy is given by the sum of energies of all visible particles: -Eghowcr = — Y^E^. 
We computed the energy spectra of final-state neutrinos in section 1231 see fig. HI (Sometimes 

decay into /i^ giving a shower accompanied by a muon: some detectors might be able 
of tagging this class of events). 

Unlike the case of /x^, shower events can be considered as fully contained at any energy. Water 
Cerenkov detectors with a high photo-multiplier coverage, such as SuperKamiokande, can identify 
e-like events most of which are due to z/g interactions and separate them from /i-like events that 
in almost all of the cases are due to CC interactions. NC and other CC interactions from 
z/g can be separated by the above topologies only on statistical basis. The energy threshold 
is lower than a GeV. Similar capabilities might be reached by a future Mton water Cerenkov 
detector. The largest and least granular planned detectors cannot distinguish from from 
hadrons: all of them are seen as showers and at the moment the energy threshold is around a 
TeV. In conclusion, it seems possible to measure the energy and the direction of the shower with 
experimental uncertainties comparable to the ones for muons. 

As in the case of fully contained /i^, we give the number of shower event for an ideal detec- 
tor with Mton-year exposure. The search in real detector requires to include the efficiencies of 
detection. In Cerenkov detectors, it is relatively difficult to see high energy z/g, so this issue is 
particularly important for showers. In the energy range relevant for the search of neutrinos from 
DM annihilation, it is possible to reach a 20% efficiency at least jlB] . 

We compute the spectrum of showers summing the three sources listed above. Indeed assuming 
that oscillations are fully known, measuring the two classes of events that we consider (yU,^ and 
showers) is enough for reconstructing the two kinds of primary neutrino fluxes produced by DM 
annihilations: and 'z^e.^t- Fig- E shows the energy spectrum of showers. By comparing it with 
the corresponding plot for fully contained muons, flg. [TJ one notices that showers have a rate 
about 2 times larger, and that retain better the features of the primary neutrino spectra (at least 
in the idealized approximation we considered). 

From the point of view of atmospheric background, and are more favorable than ^V^^ 
due to two factors: I. the flux of atmospheric drops more rapidly above E^, > 10 GeV (because 
only at low energy atmospheric /^^ decay in the atmosphere producing '■ul rather than colliding 
with the Earth). II. atmospheric are generated almost only through atmospheric oscillations, 
that at baseline 2R^ and at large energies E^ give a small P^^- ~ (40 GeY / E^Y . 

5 Neutrinos from DM annihilations in the Sun 

In this section we show the results concerning the signal from Dark Matter annihilations around 
the center of the Sun. The potential signal from the Sun is expected to be as promising as the 
Earth signal and less subject to model dependent assumptions on DM capture rates. 
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Figure 9: Neutrino fluxes generated by DM annihilations around the center of Sun. Upper row: z/g 
fluxes. Lower row: the almost equal fluxes of and v^. We plot the combination (2$j, + $p)/3. 
All spectra are significantly different from those at production point (not shown here). For the 
DM DM —7- ui? channel at Ey = m^M we plotted the survival probability rather than the energy 
spectrum, because it is not possible to plot a Dirac 6 function. 

We present our results showing the same kinds of plots previously employed in the Earth case. 
Fig.s M show the DMz/ fluxefl The main topologies of events that detectors can discriminate 
are the ones already discussed in the Earth case: fig. [TO] shows the spectra of through-going fi"^, 
fig. [11] the spectra of fully contained /i^ and fig. [12] those of showers. For brevity we will not here 
repeat the features that these two cases have in common, and we focus on their differences. 

The neutrino fluxes in fig. [n]are computed in three steps. 1) Evolution inside the Sun needs 
the formalism presented in section [3] 2) Oscillations in the space between the Sun and the Earth 
average to zero the coherencies among different neutrino mass eigenstates. The density matrix 
becomes therefore diagonal in the mass eigenstate basis. 3) Neutrinos can be detected after 
having crossed the Earth, so that we computed the functions P(z/j — z/^) and -P(z/j — z/^) (where 

These DMi/ fluxes are available at |30) . 
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Through-going fi^ from the Sun 
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Spectra of through- going fi~^ and fi~ (summed) generated by DM annihilations around 
of Sun. The plots assume the DM annihilation rate of eq. fl2^ . different DM masses 
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Spectra of fully contained fi~^ and fi (summed) generated by DM annihilations around 
of Sun. We assume the DM annihilation rate of eq. ( 124|) and a detector with Mton-year 
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Figure 12: (Idealized) energy spectra of showers generated by DM annihilations around the center 
of Sun, in a detector with Mton-year exposure and unit efficiency, assuming the DM annihilation 
rate of eq. (El). 28 



t'i = {^1, i^27 1^3} are mass eigenstates and and Ui = {ue, u^, Ur} are flavor eigenstates) taking into 
account Earth matter effects. These functions depend on neutrino energy and neutrino path. All 
plots are done assuming that neutrinos cross the center of the Earth. We also assume ^13 = 0, 
such that the actual path of the neutrino across the Earth is unimportant. Indeed for ^13 = P 
and P marginally differ from the value they achieve in the limit of averaged vacuum oscillations: 

P{ui ^ u,) ~ P(z7, -> u,) ~ \Vu\^. (23) 

Oscillations driven by Arn^^^ inside the Earth give some correction only below a few GeV. If 
instead 6*13 7^ a dedicated path-dependent computation is needed in the energy range (2 
50) GeV. 

When needed, we assume the following rate of DM annihilations inside the Sun: 

J ''^a-'^n I Earth • (2^^) 

where tse is the Sun-Earth distance, R@ is the radius of the Earth and Fannl Earth is the rate we 
assumed inside the Earth, given in eq. (12T1) . This amounts to assume that the DM annihilation 
rate in the Sun is bigger than the annihilation rate in the Earth (because the Sun is larger and 
more massive than the Earth) by a factor that precisely compensates for the larger distance from 
the source, rgs ^ -Re- SUSY models where the DM particle is a neutralino, samplings of the 
MSSM parameter space show that this compensation is a typical outcome. So our assumption of 
eq. (l24l) is again realistically optimistic. 



Backgrounds. The atmospheric neutrino background depends on the orientation of the Sun 
relative to the detector, so that we do not show it in the figures as shaded regions, and discuss 
it here. The overall magnitude of the atmospheric flux has only a 0{1) dependence on the 
zenith angle, so that in first approximation the shaded areas of fig. |5] remain similar in the 
solar case. However, atmospheric oscillations affect and ^z7^ coming from below (at least at 
energies of 0(10) GeV) and not neutrinos coming from above. The case of '^V\ is qualitatively 
important: during the day the signal of down-going solar DM'l/*- is virtually free from atmospheric 
background. Indeed, the direct production of atmospheric 'z7^ is negligible |17] and so is the flux 
from other possible astrophysics sources |18]. Unfortunately, tagging 'z7^ is a difficult task, and 
none of the proposed experiments seems able to do it. 

Furthermore there is a new background due to 'corona neutrinos': high energy neutrinos 
produced by cosmic rays interactions in the solar corona (i.e. the solar analog of atmospheric 
neutrinos) [l9]. Their flux is however of limited importance: terrestrial atmospheric neutrinos, 
restricted to the small cone of eq. (122|) centered on the Sun, remain a more significant background 
at the neutrino energies Ej^ < TeV where a solar DMz/ signal can arise. 



The effect of oscillations and interactions. As discussed in section [21 the higher density of 
the Sun mildly affects the neutrino spectra at production. Propagation effects are instead signif- 
icantly different, as illustrated in fig. [13 DMz/ from the Earth are affected only by 'atmospheric' 
oscillations: only and z/^ are significantly affected, and only below E,^ < 100 GeV. DMz/ from 
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Figure 13: Limit neutrino energy spectra at the exit from the Sun, occurring in the case of very 
energetic initial neutrinos. 

the Sun are instead affected by both 'atmospheric' and 'solar' oscillations: in the whole plausible 
energy range oscillations are averaged and all ffavors are involved. Furthermore absorption ex- 
ponentially suppresses DMu at E^, above 100 GeV. This effect is partly compensated by and 
NC regeneration, that re-inject more neutrinos below about 100 GeV. This explains the main 
features of fig. [151 

It is easier to see these effects at work looking at the DM DM — uu channel. We assume that 
the initial flux is equally distributed among ffavors, so that oscillations alone would have no effect: 
indeed in the Earth case the neutrino spectra remain a monochromatic line at E^, = moM- On the 
contrary in the Sun there is a reduced line at E^, = moM; plus a tail of regenerated neutrinos at 
El, < TTiDM- For rriDM ^ 100 GeV the line is unsuppressed and the tail contains a small number 
of neutrinos. As rriDu increases the line becomes progressively more suppressed and the fraction 
of neutrinos in the tail increases. For ^100 GeV the line disappears and all neutrinos are 
in the tail, that approaches a well defined energy spectrum. 

The same phenomenon happens for the other DM annihilation channels: the final ffux is a 
combination of 'initial' and 'regenerated' neutrinos, and the 'regenerated' contribution becomes 
dominant for niDu ^ 100 GeV. In section [5TT] we will provide an analytical understanding of the 
main features of this phenomenon. 

Table ^summarizes the effect of propagation on the total through-going /i^ rate: for ^100 
GeV oscillations give an 0{1) correction (e.g. an enhancement by a factor 3 in the case of the rf 
channel); for mDM ^ 100 GeV absorption gives a significant depletion (e.g. for ttidm = 1000 GeV 
the depletion factor is 0.1 0.01 depending on the channel). 

Table m shows the average energies of DMz/: as we now discuss solar DMz/ cannot have energies 
much above 100 GeV, a scale set by the interactions in the Sun. 

5.1 'Heavy Dark Matter' and the 'limit spectrum' 

The numerical results presented above indicate that for m^u ^ TeV ('Heavy DM') the ffuxes 
of DMi/ from the Sun approach a well defined and simple 'limit spectrum', which is essentially 
independent on the features of the initial ffuxes. In this section we describe the properties of such 
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spectrum and develop an analytic insight into them. 

There are a number of interesting cases that fall into the category of 'Heavy DM'. For in- 
stance, non-thermally produced super-massive dark matter with mass ~ 10^'^ GeV has been 
considered in [SQl E], where the 'limit spectrum' was first studied. But also thermal relics of 
strongly interacting particles yield the observed DM abundance for ttidm ~ (10-^100) TeV0 Even 
in the case of supersymmetry, moM ^ TeV is possible if the lightest sparticle is a higgsino or in 
coannihilation funnels, although this requires sparticles much heavier than what suggested by 
naturalness considerations. 

The 'limit spectrum' arises because the annihilations of very heavy DM particles produce very 
energetic neutrinos which undergo many interactions inside the Sun: the interactions wash out 
the initial features of the neutrino fluxes (that only control the overall magnitude of the final 
flux) and determine almost universal flavor ratios and energy spectra at the exit. In the Sun 
interactions are relevant at > 100 GeV and the limit spectrum is attained for primary neutrino 
energies ^ 1 TeV. In the Earth an analogous limit spectrum is attained for E^, ^ 100 TeV. 

Fig. [13] shows the outcome of a typical numerical run. The z/g : z/^ : z/,- : z/g : z/^ : z/^ flavor 
composition is determined by the effect of z/^-regeneration (which gives more *z7^ than 'z/g^^) 
and by the effect of oscillations (which equate 'z/Jj and *z7^ and generate some *z7e), essentially 
independently on the original valuesjfl The energy spectra are well approximated by exponentials 
^-E^jg ^j^]^ slopes given by ^ 100 GeV for v and S 140 GeV for v. 

In order to understand such features, it is useful to consider a simplified version of the equations 
described in section |3l which captures the main points of the full problem. Namely, under the 
assumptions that: (a) oscillations and regeneration roughly equidistribute neutrinos among the 
different flavors, so that we can replace the flavor density matrix p(-E) with a single density p{E)\ 
(b) the re-injection spectrum from NC scatterings and i/T--regeneration can be taken flat in E'^jE^ 
(although this is not an accurate description especially for z/T--regeneration); (c) the cross sections 
are proportional to the neutrino energy in the relevant energy range ((Jcc,nc oc Ey)\ the full 
equations of section [3] reduce to a single integro-differential equation with an absorption term and 
a re-injection integral: 



The spatial variable r has been rescaled here to the quantity xj S . The 'neutrino optical depth' x 
spans (0,1), where x = corresponds to the production point in the center of the Sun and x = \ 
to the exit from it. S therefore incorporates the radius of the Sun and all the numerical constants 
that appear in the cross sections. Using the explicit numbers for the CC and NC cross sections, 
we compute values for S that are in good agreement with the numerical results quoted above. 
We checked that, dropping the assumption of the flatness of the z/T--regeneration spectrum, the 
agreement actually becomes optimal. 

^E.g. technicolor models can contain stable techni-baryons (analogous to the proton) that make up Dark 
Matter [51]. TC models with a characteristic scale around a TeV were originally proposed as a natural solution 
for electroweak symmetry breaking, but now constraints from precision data push the scale at the much higher 
energies also suggested by thermal DM abundance. 

possible exception occurring if the initial flux contains only '^v\, that do not undergo CC regeneration nor 
oscillations in the Sun (due to matter suppression). We do not consider this peculiar case here. 




(25) 



31 



Eq. ( 125|) can be analytically solved: 



pix, E) = eM-Ex/S) [p(0, + ^ p(0, E') dE'\ (26) 



as one can verify either directly or passing to the variable p = e^^^'^p. The first term of eq. (126|1 
describes the initial neutrino spectrum p{x = 0,E), which suffers from an exponential suppression; 
the second term is the contribution of regeneration, proportional to the traversed portion x. 
Now, consider neutrinos with an initial energy Eq (namely p{0,E) = 6{E — Eq)): the first term 
becomes less and less relevant as neutrinos proceed to x ^ (S'/Eq; the second term (that reads 
x/(S'6{Eq — E) exp{—Ex/(S')) conversely becomes dominant over the first as x increases. More 
generally, after a path x a 'limit spectrum' of exponential shape 

p{x, E) oc exp{-Ex/<^) (27) 

is approached irrespectively of the initial spectrum p(0, E) provided that the injection spectrum 
p{0,E) is concentrated ai E ^ S jx. This is a simple non-trivial result. In all cases of 'Heavy 
DM', such conditions are verified and indeed the exponential spectra at the exit from the Sun 
[x = 1) are well visible in the outcome of the numerical computations shown in fig. 



6 Reconstructing the DM properties 

We now study how measurements of the spectra of the previously discussed classes of events can 
be used to reconstruct the DM properties: its mass and its branching ratios into the various 
annihilation channels. Even without a good energy resolution, some of the channels give different 
enough spectra that it seems possible to experimentally discriminate them. This is e.g. the case 
of uu versus bb and (to a lesser extent) versus rf . Other annihilation channels instead produce 
too similar energy spectra (e.g. W~^W~ and ZZ) so that distinguishing them seems too hard. 
This issue depends significantly on whether ttidm is in the energy range 1. 2. or 3., defined at 
page US 

Fig.s [UK (Earth case) and [TDd (Sun case) illustrate more quantitatively the discrimination 
capabilities in a specific example. We assumed a DM particle with mass mDM = 100 GeV an- 
nihilating into Tf and only, with BR(DM DM rf ) = 0.8 and BR(DM DM uu) = 0.2. 
The contours identify the regions that would be selected at 90, 99% C.L. (2 dof) with the col- 
lection of 1000 through-going muons (dashed lines) or 100 fully contained muons (continuous) or 
200 showers (dotted). No information on the rate at which this statistics of events is collected is 
assumed. An actual experiment will have an energy-dependent energy resolution, that we approx- 
imate in the following semi-realistic way: we group events into energy bins n = 1,2,3, .. . with 
energy 30 (n ± l/vT2) GeV. With a fiat probability distribution this would correspond to ranges 

^Our 'limit spectra' in fjg. llBl agree with the corresponding fig. 16 of |41j . within their uncertainties. However, in 
that and other works, such spectra arc approximated with a log-normal function, apparently with the motivation 
that the central limit theorem might play some role in determining the out-coming spectrum after many random 
interactions. We find instead that a log-normal does not fit the numerical result better than an exponential, and 
that a log-normal does not arise from the analytical argument presented above (deviations from the exponential 
form of eq. (P7)) arise mainly because the CC re-injection spectra are not flat in E'^/Ei,). 
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Figure 14: Discrimination capabilities fitting simulated data with a statistics of either 1000 
through-going muons (dashed lines) or 100 fully contained muons (continuous) or 200 showers 
(dotted) measured in energy bins of width AE = 30 GeV. The true point is indicated by a dot. 

30(n — 1/2) GeV < E < 30(n + 1/2) GeV; we assume a Gaussian probability. No energy threshold 
is assumed, and it is effectively set by the atmospheric background, that below 15 GeV dominates 
over the signal. We computed the best-fit regions that would be obtained in an experiment where 
the rate measured in each energy bin equals its theoretical prediction. In general this is not true, 
and the best-fit regions experience statistical fluctuations: we computed their average position 
(see [52] for a discussion of this point). 

The banana-shaped best-fit regions in fig. [14] arise because it is difficult to discriminate a 
harder channel from a heavier DM particle. In particular, in the Sun the spectra of neutrinos 
with energy E^ > 100 GeV are affected by absorption and regeneration that tends to wash-out 
their initial features, see section 15.11 From the point of view of the reconstruction of the DM 
properties, this wash-out is clearly an unpleasant feature. In the limit 3> TeV solar DMz/ 
approach the spectrum in fig. [13] irrespectively of the DM properties. 

Fig. [HI also illustrates that different classes of events can have comparable capabilities. Indeed 

• Through-going /i^ give the highest statistics if DMz/ have energies E^, > 50 GeV. However, 
for the purpose of spectral reconstruction, they are less powerful than contained events: 
their energy is less correlated to the energy of the scattered neutrino so that all annihilation 
channels produce similar bell-shaped energy spectra and discriminating the annihilation 
channel becomes harder. 

• Fully contained fi"^ better trace the parent neutrino spectra, but can only be observed up 
to a maximal energy determined by the size of the detector. Below this energy their rate 
is comparable to the rate of through-going yU,^. Thus, contained events are a competitive 
signal of relatively light DM particles[j 

In our figures, this can be verified by comparing fig. [HI with fig.s [7] and fig. [SI evaluating the number of 
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Figure 15: Modifications of neutrino fluxes due to propagation. The figures show the ratio of 
fluxes 'with' over 'without' the effects of neutrino propagation (oscillations, absorptions, regenera- 
tion). The lines refer to neutrinos from DM annihilations into rf (continuous line), ZZ (dotted) 
and bb (dashed), form-ou = {100, 1000} GeV (distinguishable by the corresponding maximum 
neutrino energy). 

• Showers are fully contained in all the relevant energy range (making more difficult to tag 
them) and they can efficiently trace the parent neutrino spectra, depending on how the 
detector can measure their energy. Furthermore, the showers/muons event ratio (not con- 
sidered in our fit) allows to discriminate annihilation channels that produce neutrinos with 
different flavour proportions. 



7 Conclusions 

We performed a phenomenological analysis of neutrinos of all flavors generated by annihilations 
of DM particles ('DMz/') with weak-scale mass accumulated inside the Earth or the Sun. Our 
analysis is valid for any DM candidate. Indeed, the DMz/ signal depends only on the following 
parameters: the DM mass ttidm; the DM annihilation rate, the branching ratios BR(DM DM -t- /) 

through-going muons in one year with an area A = (0.1 kni)^, as appropriate for a Mton detector. We see that, 
especially for light DM particles, the number of contained events is comparably large. 

It is useful for orientation to write the ratio of contained-to-through-going events considering fi^ continuous 
energy losses dE^/dx ~ —a (with a ^ 2 - 10""^ GeV/cm) as: 

Accent ^ a A ^ ^ / gth \ 

^through "^DM \ "^DM / 

where F is an adimensional function, and where we consider a water Cerenkov detector with threshold i?th, volume- 
to-area ratio V/A and an average efficiency of detection (e). To give a concrete example of our expectations, for a 
DM candidate of todm = 50 GeV that annihilates preferentially into taus we expect about 8 contained /i (fully or 
partially) and 16 shower events for each through-going muon event coming from the Earth (or the Sun), when we 
adopt as detector parameters i?th = 15 GeV, A = (0.1 km)^ and V = (0.1 km)'^. 
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for the various annihilation channels /. A given underlying model (e.g. supersymmetry) predicts 
these quantities: the total rates suffer a sizable astrophysical uncertainty and typically DMi^ are 
a promising DM signal. The other parameters determine the expected spectra of DMz/. We 
therefore computed the DMu signal as functions of these parameters and studied how they can 
be reconstructed from a possible future measurement of DMi^ spectra. 

We considered annihilation channels into presently known particles: uu, bb, tt, ZZ, W^W~ . 
We also considered annihilations into cc, lighter quarks and gluons: their contribution at Ej^ <^ 
is not completely negligible. Taking into account the different energy losses of primary 
particles inside the Earth and inside the Sun, in section [2] we computed the two independent 
spectra at the production point: for and for ^v\^^. These spectra are modified by propagation: 
flavor oscillations, absorption, regeneration. The necessary formalism is presented in section |3] 
and appendix O shows an example of features that simplified approaches cannot catch. Their 
combined effect is illustrated in fig. [TS] on the 'z/Jj flux for some selected cases, and amounts to a 
(9(0.1 -T- 10) correction: 

• DMz/ from the Earth are affected only by atmospheric 'z7|i ^ oscillations at energies 

< 100 GeV. 

• DMz/ from the Sun of any flavor and any energy are affected by averaged 'solar' and 'at- 
mospheric' oscillations. Furthermore, absorption suppresses neutrinos with E^ > 100 GeV, 
that are partially converted (by NC and by 'T'r regeneration) into lower energy neutrinos. 
In section 15.11 we analytically studied how and when neutrinos with energy Ey ^ 100 GeV 
approach the well-defined limit spectrum shown in fig. [131 

Our result for DMz/ of all flavors from the Earth (Sun) are shown in fig. [5] ([9]) . The comparison with 
the atmospheric background (shaded regions) is performed assuming the realistically optimistic 
annihilation rate in eq. ( |2T]) . Table Slat page 123) summarizes the average final neutrino energies. 
Table [3] at page [22] summarizes how propagation modifies the total rate of through-going /i^ 
generated by 'z/J^. 

The latter is often considered the most promising event topology for present detectors. We 
also considered other topologies of events: fully contained /i^ and showers (e^ generated by ^z/g 
and hadrons generated by all neutrinos). They have a lower rate if ^ 100 GeV but, even in 
this case, these classes of events are important because (1) their energy is more strongly correlated 
to the incoming neutrino energy. (2) there are two independent DMz/ spectra 'at origin' to be 
measured, so that at least two classes of events are necessary. 

Finally, fig. [T3| illustrates quantitatively how measuring the DMz/ energy spectra of these 
classes of events can allow us to reconstruct the basic properties of the DM particle: its mass and 
some annihilation branching ratios. 

Existing detectors and those under construction will likely not have the necessary capabilities, 
because SK is too small, and much bigger 'neutrino telescopes' are optimized for more energetic 
neutrinos (the large volume is obtained at the expense of granularity, resulting in high energy 
thresholds (~ 50 GeV) and poor energy resolution (~ ±30%)). Increasing the instrumentation 
density goes in the direction of solving this issue. If a DMz/ signal is discovered, it will be then 
interesting to tune the planned future detectors (or project a dedicated detector) to DMz/, with 
presumable energy E^ ^ IQQ GeV. 
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Note added: In the present version 5 of hep-ph/0506298 a bug in the propagated fluxes of 
antineutrinos from the Sun has been fixed, leading to corrections of the order of 10% in the 
fluxes presented in Fig. IH] and, as a consequence, in the spectra presented in Figures [TUl — )• [T^ In 
previous versions, an erroneous double counting of the prompt neutrino yield in H^-boson decays 
and a numerical bug in the implementation of the boost for top quark decays had been fixed. 
These modifications affected the W^W~ and ti channels in the fluxes at production of Figure |2] 
as well as (as a consequence) the propagated fluxes presented in Figures O —t- [12] and Tables [3] and 
lU Furthermore the values of some parameters had been updated, leading to very minor changes. 
Overall, these corrections and refinements amount to adjustments of the order of 10% to 20% 
at most in the numerical results. All physics discussions and conclusions are always unchanged. 
Updated results are available in electronic form from |30j . 
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A Neutrino spectra per annihilation event 

The neutrino differential spectrum per annihilation event is defined in the rest frame of the 
annihilating DM, since the annihilation process occurs at rest. The spectrum can be calculated 
by following analytically the decay chain of the annihilation products until a r, a quark or a 
gluon is produced. The neutrino spectrum is then obtained by using the Monte Carlo modeling 
of the quark and gluon hadronization, or r decay. We produced the = and differential 
distributions for h = T,q, c, b, gluon at various injection energies for each h {q = u,d, s stands 
for a light quark). Whenever we need the u distribution for an injection energy different from 
the produced ones, we perform an interpolation. In order to obtain the neutrino differential 
distribution in the DM rest frame we perform the necessary boosts on the MC spectra. 
For instance, let us consider neutrino production from a chain of this type: 

DMDM-^A^a->/i-^z/ . (28) 

The neutrino differential energy spectrum per annihilation event is obtained by the product of 
the branching ratios for the production of A, a and h in the decay chain, with the differential 
distribution of neutrinos produced by the hadronization of an h injected at an energy Ep^^d 
(defined in the rest frame of the a decaying particle) double boosted to the DM reference frame: 



^ = BR(DMDM ^ A) ■ BR{A a) ■ BR(a ^ h) 
dE 



ml 

~dE 



■ (29) 

boost a^Al boost A-i-DM 

The first boost transforms the spectrum from the rest frame of a (in which h is injected with 
energy -Eprod) to the rest frame of A. The second brings the distribution to the rest frame of DM. 
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Each boost is obtained by the following expression: 



'dE 



1 

2 .IE' 



'dW 



^prod 



dE' 



with E'^ = min [Ep^od, lE (1 ± f3)] . (30) 



where E denotes the energy of neutrinos, 7 and /3 are the Lorentz factors of the boost. 



B Cross sections 



The DIS NC differential cross section on an average nucleus is 



da , v-^ 2G|m 

-{uN^u'X)= ^ ^ 



dE: 



q={u,d} 
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( 1 1 E'v\ (2 2 ^'v\ 

Vn \9Lq + 9Rg^)+ Pg [9Bg + 9Lg^) 



dE: 



E?^ 



q={u,(£\ 

where Q < E^ < E^ the energy of the scattered neutrino or anti-neutrino and pu,Pd,Pu,Pd 
the fractions of nucleon momentum carried by up and down-type quarks and anti-quarks. In a 
medium that contains neutrons and protons with densities A^^^ and Np 



Pu 



Pd 



0.25iVp + 0.15iV„ 

0.25iV„ + O.lSiVp 

Np + Nn 



Pu 

Pd 



OmNp + 0.06Ar„ 
O.OSNn + OMNp 

Np + Nr, 



(31) 



The Z-couplings of quarks are 
1 2 



_ 2 _ 2 

9lu — 2 ~ ~ 3"^^' 



The DIS CC differential cross sections are 



did 



— (ued — )■ iu) = —(uid £u) 
dy dy 



1 1 2 

-2 + 3^w, 



GpS 
vr 



9Rd — 3"^W- 



da 



UiU — )■ Id) 



da 



TT 



dy dy 
where y = —t/s [0 < y < 1), The total quark CC cross sections are 



(32a) 
(32b) 



a^ued — )■ iu) = a^ued — )■ iu) = ?>a{u^u — )■ 



?,a{m Id) 



TT 



(33) 



a/S is the center-of-mass energy of the quark sub-processes. It is given by s = sx , where x is the 
fraction of the total nucleon momentum P carried by a quark, p = xP. 
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C Oscillation and absorption in constant matter 



Following [21] we described neutrino propagation by writing a differential evolution equation for 
the neutrino density matrix. Ref. [H] employed an alternative approach, merging the usual ana- 
lytical treatment of oscillations with a MC code that accounts for absorption. We here illustrate 
the non trivial interplay of oscillations and absorption, by solving our evolution equation in a 
simple semi-realistic case where analytic solutions can be obtained. Let us consider the v^/vt 
system where oscillations are due to Am^ = Aml^^ and absorption is due to matter with constant 
density. We assume maximal atmospheric mixing, an initial neutrino state with energy Ey, and 
study final neutrinos at the same energy Ey. Thus, regeneration does not contribute and the 
full system of equations (fTT]) reduces to a system of ordinary differential equations. The explicit 
solution in matrix fornJi is: 



with 



p(r) = e-"/^-'^^ ( U{r) p(0) + p(0) U\r) ) (34) 

jj ^( c-s{V^- r,)Aosc/(47r) -i s Am2Aosc/(47rE.) \ 

1^ -I s Am'Xosc/i^irE,) c + s{T,- r,)Aosc/(47r) J ^ > 

where we used the shorthands c = cos(7rr/Aosc) and s = sin(7rr/Aosc)- The factor describes 
absorption of u^^ and is approximatively given by !^ O.lAEiypGp in an iso-scalar material with 
density p. Numerically F~^ ^ 10'^ km for the typical density of the Earth mantle p ~ 2.4 ton/m^ 
and for E^, = 100 GeV. The two relevant scales that compare with the pathlength r are the 
absorption and oscillation lengths: 



Aabs = 2/(F^ + F.), Aosc = 47r/y/(AmVE.)2 - (F^ - F.)2. (36) 

Notice that absorption reduces the oscillation wavelength. When F^ — F^ is large enough Aqsc 
becomes imaginary i.e. the system becomes 'overdamped' , and oscillations disappear. While 
absorption and F^ — F^ can be neglected in terrestrial experiments (e.g., in future precise long- 
baseline experiments) their effects are relevant for neutrinos with energy E^, > 100 GeV produced 
by DM annihilation in the Sun; note that in this case oscillations into z/g are suppressed. From 
the explicit solution in eq. f p^ one finds the expressions of survival or appearance probabilities 
that resemble known things closely, e.g.. 
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